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ABSTRACT 

The stability of the flapping motion of a single blade of a helicopter rotor is 
examined using the techniques of perturbation theory. The equation of motion studied 
is linear, with periodic aerodynamic coefficients due to the forward speed of the rotor. 
So’utions are found for four cases: small and large advance ratio and small and large 
Lock number. The perturbation techniques appropriate to each case are discussed and 
illustrated in the course of the analysis. The application of perturbation techniques to 
other problems in rotor dynamics is discussed. It is concluded that perturbation theory 
is a powerful mathematical technique which should prove very useful in analyzing some 
of the problems of helicopter dynamics. 
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SUMT'IARY 


The stability of the flapping motion of i single blade of a helicopter rotor is 
examined using the techniques of perturbation theory. The equation of motion studied 
is linear, with periodic aerodynamic coefficients due to the forward speed of the rotor. 
Blade pitch feedback proportional to both flapping displacement (6g) and flapping rate 
is included . Four cases are considered: small and large advance ratio and small 
and large Lock number. The perturbation techniques appropriate to each case are 
discussed and illustrated in the course of the analysis. Analytic solutions are obtained 
for each case, with primary emphasis on the eigenvalues (that is, root loci) as indi- 
cators of the system stability and response. The feature of the equation which makes 
perturbation techniques useful is the periodicity of the aerodynamic coefficients. 

The applicability of the four cases considered is discussed; the small advance ratio 
results in particular are very useful, being valid out to an advance ratio of about 0. 5. 


The application of perturbation techniques to problems in rotor dynamics with more 


degrees of freedom or better aerodynamic models is discussed. It is concluded that 
perturbation theory is a powerful, and yet not very sophisticated, mathematical tech- 


nique which should prove very useful in analyzing some of the problems of helicopter 
dynamics. 

INTRODUCTION 

This paper considers the application of perturbation techniques to helicopter 
rotor dynamics. Perturbation theory has been well developed in recent years, but has 
not found much application to rotary wing problems. Classically helloopter engineering 






has made use of the same perturbation theories that fixed wing engineering has, for 
example lifting line theory and engineering beam theory (both require a large blade 
aspect ratio). Another classical example is actuator disk theory (a large number of 
blades is required). These theories were developed on an intuitive basis however, 
and the more rigorous mathematical techniques of perturbation theory have not yet 
found widespread use for rotary wings. The classical applications are largely for 
aerodynamic problems; the mathematics of these problems can be very complicated 
however because the equations involved are highly nonlinear partial differential equa- 
tions. The treatment of dynamic problems can be more tractable since only ordinary 
differential equations are involved. Problems with constant coefficient linear differ- 
ential equations can be solved exactly with well established methods, so for these 
problems the extra effort of perturbation theory may not be justified. On the other 
hand for problems with time varying or nonlinear differential equations the only solution 
procedure generally applicable is the numerical integration of the equations of motion. 
However, purely numerical solutions are not entirely satisfactory for obtaining an 
understanding of the physical character of the system, or for formulating general design 
rules. Furthermore, an analytic solution for the general case would be difficult to 
obtain (if possible at all) and would be so complex as to be hardly better than the 
numerical solution. The only systems that can be practicably handled analytically 
are those involving linear constant coefficient differential equations. Perturbation 
techniques are available which are methods to study time varying or nonlinear systems 
such that at each step in the analysis only linear constant coefficient equations must be 
handled. Time varying or nonlinear differential equations are characteristic features 
of helicopter dynamics and aerodynamics, primarily due to the rotation of the wing. 
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Thus ihe possibilities for the use of perturbation theory in rotary wing problems are 
very extensive. 


This paper considers the stability of the flapping motion of a single blade of a 
helicopter rotor. This is a single degree of freedom, second order system, with 
analytic aerodynamic coefficients. The governing equation is linear with time varying 
coefficients; it is given below. 

8 + l?B = y C(Mg - K r M 0 ) 8 + (Mg - K p M Q ) 8 ] (1) 


rv(i + 6 wl "* + i , " 8i " w4 ) 




M 8 = ' 


Mcosi/)^ + -Msin^ 

- /icosi/)^“+ “ Hsinip - ± (jisin0) 3 ) 
;icosi/)(~ + ~ jusin^ 


^11 . , 1 . ,.2 

-+-/ism0 + -(.usin0) 


Mg ~/isin0 + ~ (Msintf)) 2 - ~ (Msin)J>) 4 
(t + 3 M8in ^ + 4 (M sint^) 2 ) 


Regions 

(i) 

(ii) 

(iii) 

(i) 

(ii) 

(iii) 
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where the coefficients have separate definitions in three regions of the disk defined by 

Region (i) o'< p sin t/i < p 

Region (ii) - 1 < n sin 0 < 0 
Region (iii) - M < M fin lb < - 1 

This is the homogeneous equation for small perturbations of the flapping motion of 
the blade about an equilibrium state ; the derivation of this e juation may be found in 
the literature (Ref. 1). $ is the degree of freedom representing the blade flapping 
motion perturbation. The equation is nondimensionalized with the rotor speed, so 
the time variable is the azimuth angle $. u is the rotating natural frequency (non- 
dimensionalized with the rotor speed) of the flapping motion, which may be greater 
than 1. 0 for flapping hinge offset or cantilever root restraint of the blade, y is the 

4 

Lock number, defined by y = p acR /I^ (p is the air density, a the two-dimensional 

lift curve slope, c the blade chord, and R the blade radius); 1^ is the equivalent mass 

of the flapping motion, given by the integral over the span of the square of the mode 

shape of the flapping motion weighted by the mass per unit length of the blade ; for the 

rigid flapping motion of an articulated blade, the mode shape is proportional to the 

radial distance from the hinge, and so 1^ is just the moment of inertia of the blade 

about the flapping hinge. is the flap proportional feedback gain, better known as 

tan 6^; is the flap rate feedback gain. A feedback law A 6 = - K^/3 - K^/3 has been 

used (A0 is the blade pitch change due to flapping feedback control). fJ is the rotor 

advance ratio (forward velocity divided by rotor tip speed). The coefficients M«* M_, 

P o 

and are the aerodynamic forces on the blade, heuce their multiplication by y and 
their dependence on /J. The three regions for the coefficients reflect the influence of 






i 


the reverse flow region of the rotor disk. In region (i) there is normal flow over the 
entire blade span; in region (iii) reverse flow over the entire span; and in region (ii) 
normal flow outboard of r = - sin ip and reverse flow inboard. Region (iii) is 
encountered only if H > 1. The aerodynamic coefficients were obtained using a rigid 
blade motion, and should properly be changed some to handle a blade with cantilever 
root restraint. However the major effects of a cantilever root on the dynamics of the 
system are due to the change in v and y (both are increased, v to 1. 15 say and y to 
about 5/3 the Lock number based on the rigid mode inertial). Since these are free 
parameters in the analysis this formulation of the problem should give reasonable 
results for all rotors. 


This equation has been studied numerically in recent literature, primarily in the 
context of Floquet theory (Ref. 2), which must be used because the aerodynamic 
coefficients are periodic in ip if f 0. The equation will be studied in this paper using 
the techniques of perturbation theory. The mathematical techniques are not very sophis- 
ticated actually; they are rather long, especially when the higher order solutions are 


sought; and there are some tricks to be learned, but the standard ones work for most 
systems, including this one (Kefs. 3 and 4). It is not maintained that the differential 


equation studied is a true model of rotor dynamics; nonlinear aerodynamics and coupling 
with pitch and lag motions are certainly very important. The purpose of this paper is 
not to present a study of true flapping dynamics; rather it is intended to demonstrate 
what information can to obtained by the perturbation techniques, and to explore the 
methods which are most useful for rotor dynamics, so helicopter engineers will be able 
to decide whether to use these techniques with more complicated or more realistic 




systems. The dynamic problem considered here is the question of rotor flapping sta- 
bility. The stability of the motion is determined by the roots or eigenvalues of the 
system (there are two for this second order equation), and so most of the results dis- 
cussed will be concerned with the roots. The equation considered is linear; perturba- 
tion theory is used because the aerodynamic coefficients are time varying (specifically, 
periodic) for forward flight, i. e. , when M is greater than zero. A brief discussion of 
the characteristic behavior of the eigenvalues of a periodic system is given in Appendix I. 


NOMENCLATURE 

Kp Flap proportional feedback gain 

Kp Flap rate feedback gain 

Mg Aerodynamic moment due to flapping displacement 

M^ Aerodynamic moment due to flapping rate 

Mg Aerodynamic moment due to blade pitch 

0 Flap motion degree of freedom 

y Blade Lock number 

X Eigenvalue or root of the system 

M Rotor advance ratio (forward speed/rotor tip speed) 

v Rotating natural frequency of flap motion (centrifugal and structural stiffening), 
nondimenaionalized with rotor rotational speed 
Rotor azimuth angle, measured from downstream 
o( ) ’’the order of" 

0 Conjugate of a complex number 


-?- 


LHP Left-hand plane 

RHP Right-hand plane 

LHS Left-hand side 

RHS Right-hand side 

Re Real part of a complex number 

Im Imaginary part of a complex number 

ANALYSIS AND DISCUSSION 
Introduction to Perturbation Techniques 

Fundamental to the use of perturbation techniques i3 the existence of some 
parameter which is either very small or very large (how small or how large is deter- 
mined during the analysis); for the moment represent the small parameter (or the 
inverse of the large parameter) by c. In the present problem it is desired to find the 
roots of the motion, which means investigating a solution which is uniformly valid over 
long time periods. The appropriate perturbation technique is the method of multiple 
time scales . This method assumes that the behavior of the system may be investigated 
over several time scales, i. e. , 

4 ) ® c n * 

n 

The time scales are all assumed to be the same order; then for ^ * €ip the actual 
time ip must be of order € i.e. , very large compared to the basic time scale 
>Pq * ip- Next the dependent variaidos are expanded as a series in c, 

0 = 8 Q (tP 0 , 4> v 4> r ...) + c B^ipy ip v ...) + 


« • * 
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where the terms 8 ^, 3^, etc. are all assumed to be the same o.’der, and depend on all 
the time scales now. The requirement that all the 3^ be the same order for the long 
time scale behavior of the motion is crutial to obtaining the solution;, it leads, for certain 
values of the free parameters, to critical regions characterized typically by a reduction 
of the stability of the system. The details of this method will be given below in the 
context of the treatment of the rotor flapping equation. 


Often an equation of motion is such that in the limit e = 0 the order of the differ- 
ential equation is reduced. Such problems are called boundary layer problems , since 
they are characterized by narrow regions in which the solution changes greatly. The 

outer solution may be found by use of a substitution of the fo.m 

ij) 

8 = exp J p d i/) 

followed by an expansion of p as a series m f : 


P= 7 P -n + * 


p o + fp i + 


This main solution is uot valid in certain narrow transition regions or boundary layers. 
A basic part of this perturbation technique is methods to obtain solutions through tho 
transition region, so that it is possible to match one main solution to another on the 
other side of the transition region, or to boundary conditions at the base of the boundary 
layer. Again, details of the method will be given in the context of the solution of the 
flapping equation. 


For the flapping equation there are two parameters which may be used for per- 
turbation quantities: the advance ratio /i and the Lock number y. Then there are 
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four cases to be considered: small and large and small and large y. The flapping 
natural frequency V is also a parameter in the problem, but it varies little and further- 
moie always has a value at or slightly above unity (i.e. , is neither small nor large). 
Each of these four cases will be examined in turn in the following sections. 
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The Small \i Case 


2 

For the 3mall M case (to 0(M )) it is possible to ignore the reverse flow region, 
and the aerodynamic coefficients in region (i) can be used for all 0. The equation of 
motion is then (considering the case K = 0 first): 


£+(g Msin^S 


p 2 + 


/icosi/)(~ + j Msinl^ij + K p (~ + | Msin(f> + ^OXsintf) 2 )] 5 = 0 (2) 


The small parameter is the advance ratio /i; the perturbation technique to be used is 
the method of multiple time scales. The solution will be examined to 0(p 2 ). 


Hover 

For the hover case, i.e. , the limit p = 0, the equation reduces to 


which is a constant coefficient equation now. 


rhe roots are obtained from 


( 







X~£* 


(and its conjugate). 


Expansion in 


Using the method of multiple time scales, the behavior of the equation is examined 
-1 -2 

for ip of order 1, fi , ft , et< . ; that is, let 

* 0 =* 

Next expand 0 as a series in /X, with each term depending on all the time scales $ n : 

3 = 3q(0q* $ 2* + •••) + ••• 

The time derivative now becomes 


_d a_ 2 a 

dip~dip 0 + fi &p 1 + ,i &P 2 + '" 


So that the ordinary differential equation (Eq. 2) now becomes a partial differential 
equation. Furthermore, the two remaining parameters in the equation, v and y, 
are also expanded as series in Ik 


t/ = + M v 2 + ..< 


y=y 0 + M7 1 + M y 2 + 




* • • 
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(This is done because it is the characteristic of a system with periodic coefficients 
that for certain values of u and y^ there are stability degradation regions described 
by boundaries in v and y , or v and y etc.). 

Now 8, d/d</), v, and y have all been expanded as series in fJL. These expansions 

are substituted in the differential equation. It is assumed that all the coefficients in the 

expansion are of the same order; thus j/) q , ^ , etc. must all be of order 1; and 

3 Q , 8^, etc. must all be of the same order for the behavior over all the time scales 

$ n (how large is arbitrary since the equation is linear in 8, although if 8 is too large 

the equation of motion may not be valid). The equation of motion will then contain terms 
2 

of order 1 , M. M . etc. ; all the terms of like order are collected and separately equated 
to zero, to give the equation that starts the analysis at each order. 

Order 1 Results 

To order 1 the equation is 


*2 70 3 / 2 y o\ 

o 8 n +~T7~ 8. +( v' z + K_ ~ ) 8 = o 

^ Q 2 0 8 0 \ 0 P 8 / 0 


( 4 ) 


The solution of this equation is 

«o- Re Mr*2’---» eVo | 


where the root X Q is given by 


( 5 ) 


x o 2 + ^ x o + ( y o + K p?) a ° 


or 


y <> [7 


111 ) 

16 ” 0 K P 8 \16/ 


( 6 ) 
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/ ^ 

— 2 > 

. * 


and its conjugate. Since depends on all the time scales, 0^ = #q(0q> ip^, • •.), 

Eq. 4 is a partial differential equation, and only determines the behavior of 0 ^ as a 

function of ip . Thus the quantity 0 still depends on ip , ip , etc. 

0 01 1 ^ 

The order 1 equation is identical with that obtained for M = 0, i.e. , the hover 
limit, and indeed X^ is exactly the hover root (to order 1). The variation of these 
roots with y for v = 1 and several Kp is shown in Fig. 1 (the \i = 0 loci). The com- 
plex portion of the root locus is a circular arc, with center on the real axis at X = - K , 


and radius of JxP 1 + Kp 2 . The corresponding y for a point on the complex portion of 
the y locus may be obtained from the real part of X since Re X = - y/16 (no dependence 
on v or Kp). For y - 0 the locus is at X = i V (there is no effect of Kp since there 
are no aerodynamic terms if y = 0). For Kp > 0, ImX increases as y increases from 
zero; a peak in ImX is reached at y/16 = Kp where X = - Kp + i/p 2 + Kp 2 (the peak 
occurs just over the center of the circle so the frequency is given by the circle radius). 

For Kp <0, ImX decreases immediately as y increases from zero. The locus inter- 
cepts the real axis at y/16 = Kp + Kp 2 , where X = - Kp + + Kp 2 . Then as 

y the roots remain on the real axis, one going to X = - °° and the other to X = - Kp 
(the center of the circle). Thus one branch of the y locus crosses into the RHP if 
K p < 0; the crossover point occurs for y/16 * - 1^/2^; at this y the other branch 
is at X = p 2 /Kp (which is less than zero since Kp < 0). When the solution is examined 
to higher order in \k (as below), special problems occur when the frequency of the hover 
root is at or near a multiple of fc/rev. The order 1 root crosses Im X = £ for 
y/16 = K p +v£ 2 T K p 2 - 1/4. The locus will cross Im X = 1 for y/16 = Kp +7 p 2 + Kp 2 - 1. 
Since V * 1 there can be only one crossing of Im X = 1 or by the locus (except when 
V - 1, in which case the locus start at Im X * 1 for y = 0). 

-■ — — r- " 


f 

1 





The root loci for fixed y and varying P or K p would be somewhat simpler than ; 

the y loci. What is being varied is the natural frequency of the system, = vP 2 + K p ^ 
for (X complex, |X| = u)^). For the complex portions of the loci, Re X = - y/ 16 is fixed, ‘ 
so the locus would be a vertical Line in the LHP. For p 2 + K (y/8) = « the locus would 

P i 

be at Im X = ® ; for + K p (y/8) = (y/16) 2 the locus would intercept the real axis, i.e. , 

■s 

would lie at Im X = 0. For smaller P^ + K (y/8) the locus would have two branches ; 

i 

2 

on the real axis. For P + K p (y/8) = - ® the locus '*.ould be at X = ±®. The locus i 

would go through the origin, crossing into the RHP, at p + K p (y/8) = 0 (the other « 

branch would be at X = - (y/8)); since P 1 this can occur only for K p = -(8/y)P 2 < -8/y, 

i 

i.e. , for sufficiently negative K p . ; 


The preceding paragraphs have discussed the behavior of the hover root loci, 
i.e. , the fJ, = 0 roots. For /i / C, X has the same form as the hover roots, but in 
terms of y^ and P^; so it is not the entire root but rather only the order 1 part of it. 

Thus for example, if Im X^ is exactly at a multiple of i/rev for some M» this only 
implies that the hover root is near that point; the hover root, based on y = y^ + jiy^ + . . . 
and p = p q + M Pj + • • • must be a small distance (0(| n) if y^ and P^ are not zero) 
away from X^, which is based on y^ and p^ only. The order /!, ^ 2 , etc. parts of 
the roots for M / 0 will bo obtained in the analysis below. 


Order u Results 

The order M terms in the differential equation are (dropping the common factor 

of jU): 
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jL 

dti), 


8, + 


0 d 


8. + [v n +k„~ 0, 


2 1 8 yij) 1 \ 0 P 8 / 1 


2 y fy y \ 

2 ai& 1 8 o + T^^o + \T + T sin ^o)a^ fl o 


( y y y \ 

2 Vl + 'T C 0 S V‘r K P + 'r K P sin* 0 )S 0 


( 7 ) 


This equation will be regarded as an ordinary differential equation for 8^ in terms of 
</> Q . The behavior of the RHS as a function of is obtained by substituting the solution 
for 8^; this gives: 



Thi° is to be considered an ordinary differential equation for 0^ of the form 


^l d i + TW 0 8 i + ( v o 2 + K P t) 8 i 


a V 


hh + 
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where A is a complex constant (which really depends on ip , ip , etc.). The solution 
X I z 

to this equation is 


V R fii e 


*■ 0^0 _ 


1 


„ v° 

2 VT 




where 8^ e i S the homogeneous solution. The particular solution for 0^ has a 
term proportional to A^i/J^e^^O. compare this with the solution for 0^: 


* 

(constantlAj^e 


0, 


(constant)e 


’ * (constant^ 


Then 0^ will become arbitrarily large compared to 0 Q if ip Q is large enough, which 
violates the assumption that 0 Q and 8^ are of the same order for all tp^. The only 
way such a term in 0^ can be avoided as if it is required that A^ be zero. Recall 
however that the equation for 0^ is really a partial differential equation, and A^ has 
terms like d0 and 0^. Thus setting Aj = 0 gives a differential equation for 

8 in terms of ip., the solution of which carries the solution for 0. out to time 
scales of the order of n .In general, the forcing terms on the RHS of the equation 
come from the homogeneous solutions for the lower orders of the 0 expansion. It is 
the nature of the perturbation expansion (not of the particular equation being studied) 
that to each order the equation for 0 q always has the same homogeneous solution (in 


this case e ®^). Thus the equation for 0^ is being forced by its own homogeneous 
solution, which gives rise to solutions of the form $ times the homogeneous solution 
here), unless the coefficient of the homogeneous solution is set to zero. It 
is a fundamental feature of the method of multiple time scales that setting this 
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coefficient to zero gives another differential equation, which may be used to find the 

behavior of 8 .to the next time scale, 
n-l 

The method of multiple time scales, as outlined above and described in more 
detail in the literature (Ref. 3), involves then the following steps. 

a) Expand 8, d/d (!>, and all parameters as series in fX. 

b) Obtain the partial differential equation for order /i n ; write it as an ordinary 

differential equation for 8 ^ in terms of ip^; substitute the solutions obtained for 

8 , j8 , .... 8 into the RHS. 
x a n-i 

c) Find in the forcing terms the coefficient of the homogeneous solution; this 
coefficient is called the secular term . Set the secular term to zero, thereby obtaining 
a differential equation for in terms of i/j^. This is done in order that the solution 
be uniformly valid for all time . Usually it is the behavior of the solution to longer 
time scales (e.g. , 3^(0^, </^)) that is of interest, rather than the higher order cor- 
rections to the solution (e.g. , 0 (d> ), which is an order fX correction to 8 Q ). So it is 
really the differential equation resulting from the secular term that is sought. 

2 

For higher orders, OQx ) and above, the solution procedure is a bit more involved 
(there are then secular terms in the secular terms), but this is best explained by 
example. 


Returning now to the flap equation to order fX (Eq. 8), the secular term (the 

MOv 


coefficient of e 


j ) is, if X Q ± i ^ X Q : 




0 * 0 

'ox 




T 


1 


r 


or 
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where 


afl n i 


(9) 


y y 

~r (K + K J + 2UV - 77 (X A + K ) - vv 
8 x 0 P 0 1 16 1 0 P' 0 1 


X i" r 

2 Vf 


i Im X„ 


t 


The solution of this equation is 0 = Re(8 @ . . . )e^T^l); then the solution for 8 is 

01 Utt 4 


= Re ^8 


6»Re|S 02 # 2 ...)e 




+ X^ 


+ 0 Qi) 


( 10 ) 


The root is 


x * x o + " x i*-w 


+ ijv 2 + 


K £-(■£•) 
T>8 \16/ 


(U) 


f 

1 


to order M- Thus to order M the root remains the hover root; there is no effect of 
advance ratio or of the periodic coefficients. This is the case for most y and v, the 


exception being when y and v are near y Q and such that X Q ± i = X Q . 


If X Q + i =» Xq then the periodic coefficients contribute to the secular term; 
X Q + i « X Q means that y Q and are such that 

Im X 0 * i 


or 


+ K !a AT i 

P 8 \l6/ 4 
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So this case occurs when the hover root has a frequency near £/rev. (Note that 
X Q - i = is not possible because X^ has been defined to have positive imaginary 
part. ) The 0(1) root in this case is 


A 0 16 + 2 


The secular term is now 


12 V 


or 


where 


X„ + 2P v +7Kj8, 
0 0 1 8 P/ 01 


y o\ »oi / y 1 
\8 

(lo_x + — + — K )b = 
\12i X 0 12 + 6i K P/ 5 01 


_01 nR 
m>i 


Jrl v ) 5 . 

+ \192 6 K P 24 / 0 O1 


a*a r + ia i = -~+i 


2V. 




( 12 ) 


The S Q1 term arises due to the periodic coefficients. The solution to the equation 
(see Appendix II) depends on the quantity 


2 2 
D *0^ - 




( 13 ) 
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Now if D >0, has terms with time behavior like 


-(o- r ± i D)</) 1 - * ± 1 W 


and then /? Q has terms like 


e *o*„_ e -'£* ±i * (4+ '‘ D > 

P 01 


the damping is unchanged by the additional 3ecular terms, and there is an 0(H) change 

2 

in the frequency. If D < 0, 0^ has terms like 


-(0£ r iMD)^ 

3 = e ' ' 


so fly has terms like 


$ e^O-e 
P 01 


{-u ± > l ty+ i r 


there ic, an 0(M) change in the damping (both more and less stable), while the frequency 
remains fixed at i/rev. 0" » 0 must give the boundary between the two types of 

o 

behavior. Consider next the interpretation of the quantity D. Constant D implies, 
for a given y , v , and K , that a is a constant, i.e. , 

U 0 r X 

Vj - 1^(t£ - k p) - [ 1+ (t- 4K f) 2 ] 

■ ± (constant) (14) 


This equation represents a straight line of v versus y with a slope of 




yi« - k p 
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Compare this with the slope of the line represented by 



V 16 - K P 

V 16 ' "o 

2 1 
Thus the lines given by D = constant are parallel to the line given by Im X Q = — 

(these can be considered lines of v as a function of y). Furthermore, in this case 

and v give an 0(/i) perturbation from y Q and which are such that Im = ^ ; 

2 

thus it follows, since a given value of D gives two lines of v versus y^, that 
2 

D = constant represents two lines an 0(/i) distance either side of and parallel to the 

12 2 
line Im X^ = — . Two values of D are of particular interest, fi r st the boundary D = 0, 

2 

and second the maximum possible negative value of D . The latter is given by 
= 0, i.e. , 

2 „„ 

0 1 8 \16 T»/ 

This line runs through V » y 3 0 and has the same slope as the Im X » £ line there; 

11 V 6 

thus this line simply represents, to Ofyi), the Im X . ■ — line, and it is sufficient to use 

U A 

2 

the point = y » 0 for the minimum D . Return now to the solution for 0^ in terms 

2 

of D; it has the characteristics expected of a periodic system (see Appendix I). D >0 

and very large implies V and y, very large, which means a root far from Im X 

11 » 

2 

As D decreases, ReX remains at the basic value (- y/16) but there is an 0(J4) change 

2 1 
in the frequency, until at D * 0 the frequency has reached exactly Im X * — . The 

2 

boundary D * 0 occurs for nonzero values of and y^, and so the root has reached 
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1 2 
the Im X = — line while the hover root is still an 0(jU) distance away. For D < 0, 

Li 

the frequency remains fixed at i/rev while there is an 0 Qjl) change in the damping, 
both positive and negative. This type of change in the stability of the system is char- 
acteristic if periodic systems; indeed it appears here due to the terms in the secular 
equation that come from the periodic coefficients in the equation of motion; it is not 

seen in the behavior of the basic root to 0(M). There is a critical region , bounded by 
2 

D = 0, inside which the change in the damping occurs. In many problems with periodic 
coefficients, the system is unstable inside such a region; in this case however there is 
the basic (hover) damping, represented by Re X = - y/16, which is large and stable. 

The change in the damping is ±/iD, which is 0(>J) compared to the basic damping, so 
the critical region is a region of stability degradation rather than of instability. With 
this discussion as a guide, the solution of the flapping equation near Im X = — will be 

U m 

considered in more detail. 


The boundary of the critical region is given by D = 0, or 

a,, )* 

0 1 \16 P/ 8 24 v \8 P / 


(15) 


The first term on the RHS makes the line parallel to Im X = — , and the second term 
gives the width of the region; the critical region is a narrow band, of width 0(/i), about 


Im Xq = g . Outside the critical rollon there is an 0 (M) change in t^c'frequency while 

the real part of the root is unchanged from the hover value. Inside, the frequency is 

fixed at i/rev while there is an 0(M) change in the damping. The maximum stability 

2 

change occurs at the center of the critical region, where D has its maximum negative 
value, i.e. , at ■ 0. The root there is 


'r\ 






( 16 ) 
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So the 1 .aximiuu stability degradation (and enhancer if is 


MD M 
max 24 




+4K 2 

O “ V A O ‘ D Ti 


0 8 P 


which is an 0(/i) small reduction; the system remains stable because of the large hover 
damping. In general the root is given by 


X = A + a -ifi D 
Or 


where for fixed U (i.e. , p - 0) have 




Let Ay/16 = /i(y^/l6), so y » + Ay-, recall y^ is given by the requirement 

Im X Q = y must be such that Ay is 0(M) sm:>U, i.e. , y must be such that the hover 

locus is an u(M) distance from 1m X * Then 

6 


^ + 2‘ iV W 4 fe" ;< p) -»%) 1 * (fl ” 4K p) 


The critical region boundary is crossed when the quantity under the square root sign 


is zero, that is when 


M -M 




oonwr 





) 


for a fixed y (i.e. , fcr the n locus), or when 


16 ~ \ 16 / " jy .. \ 


corner 


l (m- k p) 


for a fixed M (i. ©. , for the y locus). Then the root locus is given by 




cornt. 


x " - w + 1 - 1 


corner 


) 2 - 1 


The (-i) in the last term of X becomes (±1) for u > u or Ay < Ay 

^corner r 'corner 


For Kp * 0 these expressions simplify to 


corner 16 2u 


comer 




X * * W + 2 - 1 M W V ‘fe y/Ly come/ - 1 


Furthermore, for Ji « H or Ay » Ay i.e. , far outside the critical 

corner corner 

region, the expression for the root becomes 


-t-W 


,^v 


*v 'VNlkMAK & * 
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X -W + H(m) *(«-**) 

which is just an 0(/i) expansion of the hover root when ImX is near j^/rev. Finally, 
« 2 

the solution for p^, with D >0, is 



with Of. and D given above; a similar expression may be obtained for the case 

O 

D < 0 (see Appendix II). 

These results have been used to plot the root loci for varying p and y ; the 
results obtained so far are valid for small p. Figure 1 shows typical root loci for 
varying y, with p= 0 and p = 0.1. The behavior of the hover loci (p = 0) has been 
described above. The hover loci cross ImX = ^ for y = 37. 7, 18. 9, and 5. 2 for 
Kp = 1, 0, and -1 respectively; for usual rotors then the ImX = % critical region is 
likely to be encountered only if Kp < 0. The point D on the hover locus (Kp = 1, in 
Fig. 1) is where the locus crosses ImX = i . As y increases, since this point is at 
the center of '.he critical region (Ay/16 = 0) it receives the maximum stability change, 
a; d so is pulled out to the point B. In terms of the y locus, as y increases and the 
hover locus nears Im X = ^ , the root has an 0(p) change in the frequency, pulling the 
locus toward Im X = ^ . When the locus crosses into the critical region the frequency 
has just reached i/rev, and the root locus is at the point A. For still larger y the 

c 

— — y — 1 — ’ ■ ----- — ■ — 




'I 


V 

9 . 
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frequency remains fixed while the real part of one root decreases and that of the other 
increases. When y reaches the value for which the hover root has a frequency of 
i/rev, the locus is at the center of the critical region; there the roots have their maxi- 
mum stability change (which is 0(U') so the locus is at the point B. As y increases 
more, the locus moves toward the other boundary of the critical region. The locus 
reaches that boundary at the point C, and for still larger y the frequency is no longer 
fixed at 4/rev; rather the real part of the root is the same as the hover value, while 
there is an OQx) change in the frequency which decreases in size as y increases. When 
Ay/16 is no longer 0(/i), the locus is again identical (to 0(/i)) to the hover locus. 


Figure 2 shows typical root loci for several y and varying H. The circle the 
locus starts from is the Y locus for hover QJt = 0) and the appropriate K p . The y for 
each locus may be found from ReX at M =0. since for the hover root ReX = - y/16. 

As n increases from zero, for the roots near ImX = 4 there is an 0(ji) change in the 
frequency pulling the root toward ImX = 4 . while the damping remains fixed at the hover 
value. The locus reaches ImX = 4 » i. e. , crosses the boundary of the critical region, 
when M = t* corner ' For larger fi, the frequency remains fixed at 4/rev while one locus 
moves to the left (increased stability) and the other to the right (decreases stability). 
Again, this is the characteristic behavior of roots of a system with periodic coefficients. 

2 

Order u Results 

2 

In order to find the roots to 0(/i ), it is first necessary to complete the 0(M) 
solution. After the secular terms have been removed, the equation for fi becomes, 
for ImX^ / 4: 





T 
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= -~Q~ &QI « X o + 2K p^ Sin ^0 + COS e ^°^° + conjugate (24) 

The solution of this is 

^ = Re p n (0 X ) e X 0^0 + /3 qi (A 1 sin 0 Q + A g cos 0 Q ) e X °^ 0 j (25) 

where 

X 0 + 2K p • 2i Im X 0 1 * < X 0 + 2K P> 2i Im X 0 

A - 2 ’ ^2 = 9 

1 - (2 Im X Q ) z 1 - (2 Im X Q f 

or recalling that /J^ = jS^ e X l^l, the solution so far is 

V' Re (W 6 * 1 * 1 eX °*°j 

^H#l) e^ 0 + ~ P Q2 (0 2 ) (A x sin 0 Q + A g cos 0 Q ) e X l^l e X °^°] ( 26 ) 

2 

The order /i terms in the equation of motion are (dropping the common factor 
of IX 2 ): 




,2 y 


vf *: v( , '* 2+,t p ifo 


* 2 «1 + 'f^[' S 1 *(t * T sin *o) 5*;' s i 

y y y \ 
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\ 8 6 0 / 3^ 0 0 
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♦ l^i 2 + a v* + 7" 008 ♦o + f 3in 2 * + K v\f+f sin *0 + f (3 ‘" V )J *0 


h !o\ 

*« 

h 8 1 

a* x 

y 

+i 

sin 0 Q 

. y s 

(M 

oa° 

/O 


(f + f sto O x o + 2 Vi + f coa *o +K p(‘r + f sta *o) ? u eX °*° 


(\ 2 + (lf + y f a,n *o) \ + (j - + '■» S ‘“*o) X < 


y n y. 


0 8 / 00, 


2 r i y 0 

+ v. + 2ia u + — cos w n '+ — - sin 2 </) 

1 0 2 6 M) 8 0 

+ K p('T + "»• 8in *0 * f (8in V 2 )) s 02 eXl * 1 ^ 

y r y y 

* t S 02 [( 2 V'T + f 8 ln *o )(< X 0 A 1 - A 2 > 8ln *o +(X 0 A 2 +A l )cos *o) 

*(t x i ♦ a Vi + ^ 008 *o + k p(it + if 81 n *o)) <A i aln V A 2 008 V 
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Considering now Im / 1 (already have assumed Im X^ ^ |), the secular term is 


^-Vu)=^(-V 




2 i Im X. 


*02 <M (28) 


where X^ is the 0(p ) term in an expansion of the hover root, i. e. , 
x » * " x i * * - w ■ * o <" 3) 


Regarding this is an ordinary differential equation for fi in terms of $ , the RHS is a 

constant times the homogeneous solution. In order that the solution be uniformly valid, 

i.e. , that 0^ be no more singular than 0 , the secular term of this equation must 

also be set to zero. T^he result is a differential equation for j3 in terms of if) , which 

2 

will give the root to 0 (/x ): 


r X 2 + ( K P T + (t) -T 1 + A lV< 2 1 Im vK * 0 


Thus the root is 


x - x o + " X 1 + " 2 ( x 2 - ( K P T + (t) -T- 1 ♦ k p(t) a i) / < 2 1 Ira x 


, J ( y °) 2 **** '' IK 

Vor \ U 4 h » t (.-(^ v 1 P 16 Im X < 


To order jx this result is 


^ +i ^ 2 + < 1+ M 2 ) J 




?/. ..2 8 ^-8 K P + 4 V 


, 1 -<* 5 


(Z-f.J.Zv + i 

\16/ 8 F 4 





> 


Thus to order ju the roots for most y Q and v 1 (away from Im X = k or 1 that is) are 

2 

just the hover roots with an 0(/i ) change in the frequency. There are two effects of 
M; the first corrects the term y/8 Kp in the hover frequency to properly account for 
the average of K p Mg, i.e. , multiplies this term by the factor (1 + fi ). The second 
effect, that in the last term of the frequency, is entirely due to the periodic aerodynamic 
coefficients; this is the first effect of the periodic coefficients seen in the analysis, 
except for the critical regions near Im X = Typical root loci for varying jU , con- 
structed from Eq. 30, are shown in Fig. 2. These are the loci that are not near 
Im X = k or 1; the frequency change is small even at ju = 0. 5. Equation 30 may also 
be used for the branches of the root loci on the real axis when the quantity under the 
square root sign is negative (i. e. , for y large enough). There are two real roots then, 
the (+i) in the frequency becoming (±1). A point on the locus of special interest is where 
one branch of the locus on the real axis crosses into the RHP, i. e. , becomes unstable 
(see Fig. 1). The criterion for this divergence boundary is that X = 0, or 


p 2 + (1 + /i 2 ) K = - (^-) fi 2 — 
' H ' 8 P \16/ H 9 


2 y 2 

* -8 K P + 4K P 

(xi* V_z K+i 

\16 ) 81> 4 


To order j Li this becomes (since y/8 K p must be an 0 (|i Z ) distance from -p 2 ) 


/ \2 2 
(lS. V_ 

+ (l + u 2 ) ^ k =-p 2 ju 2 ~^i — L 

1 M ' 8 P M 9 , v 2 . 


m + i 

\ 16/ 4 





i^jn * 


(’.-m 2 )|k p . 


2 L 2 16 1 

-1/ (l+M g- 


The effect of M on the RHS (due to the periodic coefficients) dominates that on the LHS 
(due to the average of K p M^) for all values of y and v. Thus the critical value of 
negative Kp, beyond which the locus lies in the LHP, is actually increased by increas- 
ing fJL. The criterion from the hover case is conservative then; this is the opposite of 
the conclusion that would have been reached from a consideration of the averaged 
coefficients only. 


If X Q + 2 i = X Q then the periodic coefficients contribute to the secular term; 
Xq + 2 i = Xq means that v and are such that 

Im X = 1 


V 0+ K ? 


/V| 

8 V 16 / 


Thus there will be an 0(j it ) critical region where the hover root is near Im X = 1 The 
0(1) root in this case is 


V-Te + I 


t 
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The secular term for the equation for j3 Q2 is then 



The first term on the RHS is the same as the RHS of Eq. 28 (with Im = 1). Unless 

X is real, so that X = X , the secular equation for 0 will therefore be exactly the 
111 0 * 

same as above (Eq. 29). Recalling that 

*1 , y i( y o \ 

*1 " Um X 0 16 Im X 0 


Im X » 0 requires 


2VV 3(^. k ).« 

o i 8 Vie p/ 


This is recognized as a line of constant Im X; since it goes through 0, it is 

Just the Im X » 1 line to 0(M). and it is therefore sufficient to consider only the case 
v » m o (so Xj » 0 also). Then the secular term in Eq. 32 is 


#. 


-KK +! V2 +K P £)♦(£) u( V 0 K sWK M 1 


t. 
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The solution depends on 



The behavior of the solution near the critical region is similar to that near = 

2 2 
The boundary is given by D =0, which gives a narrow band, of width 0(/i ) here (as 

opposed to 0(/z) for the Im X Q = ^ case), about Im X = 1 . Outside the critical region 

2 2 
(D > 0), the damping is the same as the hover root and there is an 0(/i ) change in the 

frequency; at the boundary of the critical region the frequency reaches l/rev. Inside 

2 . 2 
the critical region (D < 0) the frequency is fixed at 1/rev while there is an 0(^4 ) 

change in the damping, with one root becoming less stable and the other more. 


The boundary of the Critical region (D 2 = 0) is 



This is a line parallel to Im X * 1. The terms on the LHS give a line parallel to 

2 

Im Xq = 1; the first two terms on the RHS giv* the correction required due to the Ofyi ) 

change in the frequency found above (Eq. 31), so the line is parallel to Im X - 1 for the 
2 

basic root to 0(M ). The last term on the RHS gives the width of the critical region; so 

2 

the critical region is an 0 (M ) band around Im X » 1. 





r 
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i.e. , 


2 

The maximum stability change occurs for the maximum negative value of D , 
at 





K p + 4K p 2 
12 


+ K 


P 16 


( 36 ) 


2 

which is just where Im X = 1 for the basic root (the hover root plus the 0(M ) correction 
to the frequency). At this point the root is 


X = i 



D 

max 


So the maximum stability degradation (and enhancement) is 



2 

which is an Ofyt small change; the system remains stable because of the large hover 
damping. In general the root is given by 


Vw -‘ D '‘ 2 


where for fixed v (i.e. , 9 0) we have 





(39) 




-34- 


Let (Ay/16) = M {>' /16), so y *y + Ay; recall / is given by the requirement 

u v U 

2 

Im X^ = 1; y must be such that Ay/16 is 0(m ) small, i.e. , y must be such that the 

2 . 

hover locus is an 0 (fi ) distance from Im X = 1. Then 


where 


# ‘ 2 D > 2 = (l?(u - K p) - " 2 C l) 2 - <" 2 °2 


y 2 -“K +4K 2 y 

— + K plf 


C = _(10) § — E 

1 V 6 / 12 


v!^t(iT 2K p)] + [ K P-^-2K p ) 2 ] 2 


2 2 

The critical region boundary is crossed when <JX D) =0, that is when 


2 2 

M =M 


16 2 16 P/ 


2 2 


corner *C ^1 * ^2 


(40) 


for the M locus, or when 


Ay 

16 


fe) 


M (C, ±C 0 ) 


corner 


“ - K 
16 P 


1 *°2* /Ay\ /Ay\ 

.. 'WrH 


(41) 


for the y locus. Then the root locus is given by 


X - 16‘ + 1 - 1 (w)(w ' “p) ' mVp, 2 H1 - t»V 




(42) 
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or 

X = - ^ + i - i n 2 y Cl 2 - C 2 2 ^Ay/Ayj - 1 ) (Ay/Ay,, - lj (43) 


The (-i) in the last term in X becomes (±1) inside the critical region. These results 
have been used to plot in Figs. 1 and 2 typical root loci near Im X = 1 for varying y 
and /i. 

Flap Rate Feedback 

The use of flap rate feedback, K = 0, does not change the behavior of the 

R 

solution qualitatively. The hover root becomes 



and there are cvitical regions about Im X^ * i and 1 again. The critical region boundaries 

and stability degradation depend on now. It is 

root to be stable, but K > 0 will be the usual case anyway. 

R 

y - u Plane 

The results of the small u analysis may be usod to plot lines of constant Re X and 
Im X on the y - M plane. Typical results are shown in Figs. 3, 4, and 5 for V » 1. 0 
and Kp ■ 0, 0. 1, and -0. 1 respectively. The critical regions appear in the y-M plane 
as regions in which Im X is constant &/rev or l/rev)*, they are indicated in the figures 
by the circled values of Im X (the region where Im X ■ 0 is where there are two real 
roots, not a critical region). These figures are interpreted as follows. A horizontal 
line is a line of constant y, and so as ji varies it gives the corresponding value of Re X 


necessary that K R > -1 for the hover 
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and Im X as a jU root locus does. Similarly a vertical line is a constant fJi line, and 
so gives X as a function of y just as a y root locus does. For example, consider a 
horizontal line in Fig. 3 (Kp = 0) v/ith y = 8, i. e. , y/16 = 0.5. As H increases, the 
line remains parallel to the Re X = constant linos so Re X remains fixed at the hover 
value. The Im X = k region comes closer to the horizontal line as fi increases, which 
means that Im X moves toward i/rev. Eventually the constant y line crosses into the 
Im X - region; then Im X is fixed at i/rev while for each point in the region there are 
two values of Re X, giving the damping for the two branches (one more and one less 
stable than the hover root). This behavior is just that seen already in the /J loci 
(Fig. 2). Figures 3, 4, and 5 may be compared with similar ones in Ref. 2, which 


were constructed from numerical calculations; on the basis of this comparison, the 

• 2 

0 (n ) analytic results are quite accurate up to \i = 0.5 or so. There is some discrep- 

1 

ancy between the results for the Im X * X region however, particularly with Kp = -0. 1, 
although the change in scale (Ref. 2, shows results out to » 2. 5) exaggerate the differ- 


ence. For v exactly 1 the analytic results indicate no Im X * 1 critical region if 
Kp < 0 (for v Q * 1, Eq. 35 shows that the critical region at y^/16 » 0 has zero width 



* 

1 


\ 


unless Kp = 0); but only a slightly larger v (for example v ■ 1.01) is necessary to get 
a sizable critical region with Kp * -0. 1 (see Fig. 5). The analytic results show the 
v = 1 case is a very sensitive one for small y, and it is unlikely that a numerical calcula- 



tion would treat the case accurately. Of course an actual rotor will always have V at 
least slightly greater than 1, so the numerical calculations would be reliable then; 
furthermore, the discrepancy may also be an indication that lor very small y the 
analytic results are not valid out to as large a Ji as they are for more reasonable y. 
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In any case this discussion illustrates the kinds of problems that may be hidden in a 
purely numerical solution; 'they can only be found and studied by analyti procedures 
(which at least tell where to look for problems). 

Reducti i to Mathieu's Equation 

The equation for rotor flapping stability may be studied by converting it to 
Mathieu's equation; Mathieu's equation is an equation of the form 

A 

— jr + (a - 2 q cos 2 z) y » 0 
dz 2 


It is the classic example of a differential equation with periodic coefficients, and the 
functions satisfying xt (the purely periodic ones are called Mathieu functions) have been 
well studied and documented. To transform the flapping equation to Mathieu's equation 
it is first necessary to remove the /9 term; this is done by the substitution 

p ® ye 

This is equivalent to separating out the hover damping from the solution. Then the 

classic instability regions of Mathieu's equation (for certain values of a versus q) are 

just the critical regions without the large (stabilizing) hover damping. To get the flap- 

2 

ping equation in tbe required form it is also necessary however to neglect all 0 (p ) 
terms; when this is done one obtains (with Kp * ■ 0) 

\ 2 ] 


{•Mw). 
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and 


2 z 


= $ + tan 



Thus using Mathieu's equation implies only an 0(|U) analysis. There are of course 
classical techniques for handling the more general equation with periodic coefficients, 
of the form 


^ + f(2z)y = 0 
dz 2 




which is called Hill's equation; f(t) is a general periodic function with period T = 2 w. 

These techniques could be used to study the flapping equation, to all orders in fl. 

However these general techniques are all unsatisfactory in that they tend to obscure the 
physics of the system being studied, both because some transformation is necessary to 
arrive at the required form of the equation, and because using standard solutions or 
formal calculation techniques means the great amount of information gained in the process 
of deriving the solution is lost. Furthermore, the classic treatments have only con- 
sidered a single degree of freedom system, so that they are not immediately applicable 
to more general problems. 


The Small y Case 


£ 


Consider for the small y case the flapping equation with both proportional and 
rate feedback. Now /i is arbitrary, so the general equation is considered, of the form 


W. 


3L 


$+ V jS = y C(M£ - K r M 0 ) b + (M 0 - K p Mg) j8] 


(45) 

2> 






v|| 


where the aerodynamic coefficients are functions of M and ijj (periodic in ip), and 
include the effect of the reverse flow region (see Eq. 1). The small parameter is the 
Lock number y. The perturbation technique to be used is the method of multiple time 
scales. The solution will be examined primarily to 0(y). 

Zero Lock Number 

In the absence of aerodynamic forces, i.e. , the limit y = 0, the solution is 

0 = Re (fij, e 11 ) 


so the roots are 


X = iv 


(46) 


and its conjugate; i.e. , the solution is an undamped oscillation at the rotating natural 
frequency of the flap motion. With no aerodynamics, there is of course no effect of /!. 
This result for the root agrees with the low fi result for y = 0 (see Fig. 1). 


Expansion in y 

Using the method of multiple times scales, write 


* 


0 

1 


=y<l> 


2 , 
=y <l> 
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and then expand 0 and (p/<pij) as series in y: 


0 = 0 Q • • . ) + 7 # Q , • • • ) + 


d _ 5 a 2 9 

d >l> ~ Z$ 0 +y W 1 +y W 2 + 


Also expand the free parameter v as a series: 


v = V Q + y U 1 + . 


Order 1 Results 


To order 1 the equation is 


3 +u 2 8 =0 
* , 2 p 0 0 P 0 

% 


The solution of this equation is 


V Re 


W * 2 - : 


‘Vol 


So the roots are 


\ miv o 


and its conjugate. This solution is of course just the y = 0 limit. 


Order y Results 


The order y terms in the differential equation give (dropping the common factor 


of y): 


t 
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+ conjugate 


(49) 


Now expand the aerodynamic coefficients as complex Fourier series: 


V2 

n=-°° 


where 


217 

e M a 0 * 


and similarly for M. and M.. 

P o 


If + n ^ i. e. , ^ n/2 for any integer n (i. e. , Im not equal to a 

multiple of i/rev), then the secular term has no contributions from any harmonics of 
the aerodynamic coefficients except the zeroth harmonic (the average over the azimuth); 
setting the secular term to zero gives 


2f + [-‘V • Ve) + 3 tH- K p“e)] 


B = o 
p 0l 


(50) 
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The solution of this equation is 


8 = 8 (ib . . . ) el 

P 01 02 W 2 ’ 




(51) 


and the solution for 8 is 


0-Re « 02 e 


♦ r [« v * K - V4) - 4 K - K P M J 




>+ 0(y) 


(52) 


The roots are 


X = i V Q + y 


lv i + i ( M r K R M e)-itK-Kp M e) 


and its conjugate. Before proceeding further, it is noted that for the particular aero- 
dynamic coefficients considered here many of the harmonics are zero. Actual calcula- 
tion of the harmonics or symmetry arguments can demonstrate that 


lc 2s 3c 
M^=M£ =M£ =... =0 


0 Is 2c 

m” = Mg = Mg = . . . = 0 


..lc --2s .,3c „ 

M e = M e = M e = ••• = 0 


(53) 


(the superscripts relate to a cosine/sine Fourier series). Using the fact that = 0, 
and also that to 0 (y) v = + yv , the root becomes 


x = 2( M r K R M e) + - ( 1 + ^ K P M e) 




-- 


(54) 
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Flap rate feedback only effects (to 0(y)) the damping; the real part of X may be written 


Re X = £ Ui 

2 0 


1+ V 


/ ^1 

U»)J 


The ratio -M^/Mg determines the relative effect of K R ; this ratio is shown in Fig. 6 

as a function of /A. This parameter is a positive number, which varies little with [i 

(from 1 at M = 0 to £ at \i = <=, with most of the change below \l = 1). The negative 

of this ratio gives a critical K , since it is necessary that K > K = M^/M° 

R R crit * 6 

in order that the system be stable for small y. K R > insures stability for all jU 

(to 0(y) and with v q t n/2); for very small M the criterion is K R > -1, which agrees 

with the result from the hover root (Eq. 44). For K = 0 the root is 

R 


V 1 

*r- 

> 


x= ~u 




(55) 


The aerodynamic coefficients (-8M^) and (8M^) are always positive; they have the 
value 1 for - 0 and are asymptotic to 8/37T M and 16/37T U respectively for large \i ; 

o 

these coefficients are shown in Fig. 7. To order fl the root is 

2 
V 


X --W + i, { 1 + £-2 < l 

•* V 


+ M 2 ) 


which agrees with the 0(y) expansion of the small H results for X (to 0(/i ), but with 
or 1). 


The root loci for varying y and varying H are shown in Fig. 8 for K p » 1, 

K r = 0, and v ■ 1; the locus for K p * -1 is obtained by reflecting this locus about the 

— r— • , ~ 
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Im X = v line. The locus for K = 0 is difficult to plot since Im X = v for all y and /i 

r 

(to order y); it may be visualized by projecting the K p = 1 loci onto the line Im X = v. 
These loci should be compared with the small y portions of the curves in Fig. 1, which 
are for small \l. In Fig. 8, the y locus for a given fx starts out at X = iv always, and 
is a straight line with slope 

Mm A K P 1 
aReX “" * 

which varies from -K p /t' to -2K p /v for \x from 0 to °°. The step size on the y locus, 
for a unit change in y/16, is 

+ (K p /V) 2 (^8M°) , 

which varies from 

A + (K p M 2 to f/l + f <K p M 2 t0 M M A + 4 (K p M 2 




for fi from 0 to 1 to * respectively. The M locus for a given y starts out vertically 

from the M = 0 line, and is asymptotic to the M — 00 line, with the step size on the locus 

for a unit change in ju increasing as H increases. For reasonable ju the locus does 

not vary much from the small M results. An 0(y) analysis can only obtain the slope of 

the y locus at y » 0, so the locus is a straight line, as found above; to find the curva- 

2 2 

ture effect it is necessary to go to order y . The significance of the curvature (0 (y )) 


may be judged from a comparison of the y loci of Fig. 8 (all M> small y) and Fig. 1 f 

(small M» all y); on this basis the small y results should be limited to y/16 =0.2 | 

\ 



A WesaGiZ 


or less. On the basis of neglect of the curvature effects alone the results might be 
accepted to higher y, but the 0(y) results will also be limited by the effects of the 
critical region, which will be examined next. 


The discussion of the 0(y) analysis has so far only been concerned with the basic 
roots, meaning the roots away from the influence of a critical region. In this problem, 
the criterion for being away from a critical region is that / n/2 for any integer n; 
this may be written v / (n/2) + 0(y), i. e. , the rotating natural frequency may not be a 
distance of order 0(y) from (n/2)/rev. Since V is almost always just slightly above 
l/rev (V = 1. ?. would be very large for a rotor; it would require very stiff blades and 
thus also mean high blade loads) this criterion is seldom fulfilled, and the critical regions 
may be expected to dominate the root loci behavior for small y. Furthermore, if Kp 
is large enough positive or negative, the basic locus will also cross Im X = 3/2 or \ 
for y/16 still small (see Fig. 8), so these critical regions may affect the loci even if 
the Im X = 1 region does not. 


If - p + n = v for some integer n, then the higher harmonics of the aerodynamic 

coefficients contribute to the secular equation, and there arise critical regions with 

behavior of the root loci similar to that encountered already in the small M case. The 

criterion - v Q + n * means U Q * n/2, i. e. , V * V Q + y v + , . . is 0(y) from a 

multiple of i/rev ; the only cases likely to be encounteref for rotors are v Q = fc, 1, and 

3/2 (n * 1, 2, and 3). Only the case K » 0 will be ..j Hide red for now, and use will 

R 

also be made of the fact that » 0. Setting the secular term in Eq. 49 to zero gives 
then 





T 
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■[K) + 4K 


p p 


p M e 

j *01 

■ M e)" 

'•d 



where recall that M 11 is the nth harmonic in the complex Fourier series expansion of 
the aerodynamic coefficient, and n is here given by n = 2i/ . For this equation 

< 57 ' 

Recall that for the aerodynamic coefficients considered here, either the cosine or sine 
term for each harmonic is zero (Eq. 53) , so this parameter may be written 


° 2 - { v i + M e) - (K - ‘ 'tfL 2 + K P 2 l M ef ) (58) 

For n odd, m” - i and m” = - ~ M° S ; for n even, 

* ” 2 (^ 8 + ^0 = ^ M 0 C - The boundary of the critical region 

2 2 

is given by D = 0 ; outside the region (D > 0 ) the real part of X remains at the basic 
root value (y/2 M^) while there is an 0 O') change in the frequency; inside (D < 0) 
there is an 0 (y) change (both positive and negative) in the real part of X while the fre- 
quency remains fixed at (n/ 2 )/rev. Constant O means ( 21 /^ + K p Mg) » constant, or 

dp " K P w° 




Compare this with the slope of the basic root (Eq. 55) when Im X = constant: 


so 


Im X = v 4 


Z K P 0 

t* ~ M = constant 
^ 1/ o 



rh 

2V„ 


M 


0 

0 


So the critical region is a narrow band, of width 0(y), around Im \ = (n/2)/rev. The 
2 

boundary (D = 0) is 


o M d| 2 + K 


M, 


(59) 


The maximum stability change occurs at the center of the critical region: 


v 


1 




which is where the basic root would cross Im X = - n/2; there the root is 


x " i V2 m S1 ±, ' d 


max 


‘V 


2( M fl ± ^4 M S- U 'oS| 2 + K p 2 | M 9f) 


(61) 


Inside the critical region the frequency is fixed at = n/2 whhe there is an 0 (y) change 
in the damping. The damping of the basic root is itself 0(y) however, so in contrast to 
the small \t case, the critical region can here lead to actual instability, not just stability 
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degradation. In general the roct is given by 


\ Z--° 

X = 1 + f- 


o ■ 2 M e +lyD 


( 62 ) 


Let Av = yv , so v = v Q + Av = n/2 + Av and v must be such that Av is 0(y) small, 
i.e. , an 0(y) distance from \'n/2)/rev. Then 


X = 


'V 




(63) 


Far outside the critical region this result approaches the basic root (Eq. 55); however 
for the small y case it is difficult to get far away from all critical regions, since 
|Ay| = \v - n/2| is at most — which is not very small. The critical region is crossed 
when yD = 0, i.e. , when 


y=y 


2v q Av 


corner 


"W* 




+ K„ 


ni 


= y i* y 2 


(64) 


M. 


So the y locus is given by 


X = i v Q + 1 + i Av 7(1 -y/y^l -y/V 2 ) 


(65) 


Inside the critical region the (+ i) in the last term in X becomes (±1). The M locus is 
best found from Eq. 63 directly, since the harmonics of the aerodynamic coefficients 
are rather complex functions of 





The H loci show the behavior characteristic of periodic systems, and familiar 
from the discussion of the loci for small /i. For small ji, Re X is fixed at the basic 
root value while the frequency moves toward (n/2)/rev. For v near 1 and small y, 
the locus moves toward Im X = 1 for small Kp, toward Im X = 3/2 for Kp = 1 or so, 
and toward Im X = £ for K p = 1 or so. For some the locus crosses the boundary of 
the critical region; at this point the frequency has reached (n/2)/rev. For larger M 
Im X remains at n/2 while the effect of the critical region is to decrease the stability 
of one root and increase that of the other. The maximum stability change occurs at the 
center of the critical region. The center is reached when (see Eq. 60) 

<66) 


The aerodynamic coefficient (8M^) is a positive number greater than 1 and monotonically 
increasing with ji; it is shown in Fig. 7. If / 0, this criterion will always be satisfied 
for some y/ 16 , so the y locus always reaches the center (just as the y locus for small 
/i always goes through the center — see Fig. 1; however in this case the value of y/16 
required may be outside the range of validity of the solution if Kp is too small). For 
the locus, again if Kp / 0, the center will be reached for some M provided At/ has 
the same sign as -Kp, and |Ai/| > |Kp|/(n/2) y/16. However, if Kp = 0 the center is 
never reached for either the y or n loci unless At/ ■ 0, i.e. , v is exactly n/2; in 
that case the locus is always at the center, since the frequency of the basic locus is fixed 
at v ~ n/2 then. This discussion has just been a more quantitative examination of the 
criterion that the center of the critical region is reached when the basic root would have 
crossed Im X = n/2; it is illustrated graphically in Fig. 8 for Kp » 1 and n/2 3 3/2. 






/ . 
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As such, the criterion is limited by the fact only an 0(y) analysis has been used; the 

curvature of the root loci due to 0(y^) effects can be quite important, particularly for 

small K p . For example, for K p = 0 the basic locus has a slope of zero to 0(y), i.e. , 

2 

Im X is equal to v for all y and H, and in general for sou 11 Kp the 0(y ) change in 

2 

the frequency will be more important than the 0(y) change. An 0(y ) analysis would 
change significantly the conclusions about whether the center of the critical region is 
reached under certain conditions; for example, Fig. 1 indicates that with v slightly 
greater than 1 the y locus (for small fi) would never cross Im X = 3/2 for K p = 1 
while it would always cross Im X ■ 1 for K p = 0, just the opposite of the conclusions 
indicated by the 0(y) results (Fig. 8). In any case, since the maximum stability 
change occurs at the center of the critical region, it is useful to examine it as a worst 
possible case, which may perhaps be appr oached but never reached for certain values 
of v. 


Returning now to the behavior of the M locus, the maximum change in the damp- 
ing from the value of the basic root is (Eq. 61) 



. 0 

The contribution from the basic root damping, (y/2) M^, is always negative (Fig. 7); 
as for the small M case it is 0(y), but here that means the basic damping is small. In 
fact it is the same order as the critical region contribution, so the destabilized root 
may be actually unstable, rather than just a small perturbation from the basic damping 
as for the small H case. The behavior of the locus depends on the relative effects of 
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Mi and the nth harmonics of the aerodynamic coefficients under the square root in 

P 

Eq. 67. For most cases the critical region effect dominates, so that as p increases 
it eventually reaches a critical value, at ’ hich point (for the case of the maximum sta- 
bility change) one root crosses int^ RHP, i. e. , becomes unstable. From Eq. 67, 

2 

it follows that increasing (Kp) (Kp either positive or negative) always increases the 
effect of the critical regions, which means decreasing the critical p for which the root 
becomes unstable. Thus the critical p is a function of |Kp| , for each of the critical 
regions; this function may be found from Eq. 67 by setting Re X = 0 /the requirement 
for crossing the Im X axis). Since the aerodynamic coefficients are rather complex 
functions of p, it is more convenient to find the critical |Kp| as a function of li: 



This may be regarded as a maximum K_ for a given p; for large K. 


the locus is 


in the RHP at that p. These boundaries of |Kp| versus p are shown in Fig. 9 for 
Im X near £, 1, and 3/2. With the exception of roots near Im X = £ (which requires 
Kp < 0 since v is near 1) with j- above C. 5 or so, Fig. 9 shows the criterion on Kj 
is not very stringent; a value for |Kp| of 2.0 for example is quite large, correspond- 
ing to |fi 3 | * 63.4 degrees. Figure 9 shows also that near Im X * 1 the roots are always 
stable, regardless of p, if |K p | < «/2; the loci may be expected to be near Im X = 1 
for zero or small | | . In terms of the p locus, this means that as p increases the 
locus does not cross into the RHP. Just after the locus crosses the critical region 
boundary, the effect of the critical region is seen and one branch moves to the right and 




~7 


s . 



k 
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the other to the left (as do the loci in Fig. 2). As fJ. increases further however, the 
damping of the basic root (which is always stable, and increases with /i) eventually 
dominates the effect of the critical region, and the root which was becoming less stable 
turns around before reaching the RHP. So for larger fJ, both branches of the locus will 
be moving to Iiie ieif, i. e. , becoming more stable as H increases. For the root loci 
near Im = \ or 3/2, the effect of the critical region remains dominant, and so one root 
eventually crosses into the RHP as ju is increased. The critical /i is considerably 
lower for Im = | than for Im X = 3/2. This points out an undesirable feature of negative 
pitch flap coupling, Kp < 0: not so much that it reduces the critical fl, bu* rather that 


it moves the basic root nearer to Im X = 


Flap Rate Feedback 

The use of flap rate feedback, K / 0, results in no qualitative changes in the 

R 

behavior of the loci. K is however a useful design parameter; it may be used for 

R 


example to raise the critical /i or 


K 


Evaluation of the Order y Results 

Numerical calculations were made of the fi root loci for moderate and small 

values of y. On the basis of a comparison of the numerical and analytic results, it 

is concluded that the small y analysis to order 0(y) is useful only for truly small y, 

e.g. , y = 2 or 3 (y/16 = 0. 2 or so). Problems are encountered with both the basic 

roots and the effects of the critical region. The basic root to order y neglects the 

curvature of the y locus, which is especially important for zero or small Kp, since 

2 

then the change of Im X for small y is due more to 0(y ) terms than to the 0 (y) term. 
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a 


For example, a locus based on the 0(y) analysis would start out (Jl = 0) at the 

wrong point, the error being the difference between the circle giving the exact y locus 

at = 0 and a line tangent to the circle at y = 0 (see Figs. 2 and 8). The damping of 

the basic root is 0(y) always, no matter what order the analysis is carried to; e.g. , 

2 

the 0(p ) results give Re X = - y/16 for all y. But while the basic damping is 0(y), 

2 

the contribution to the damping due to the critical region will have terms that are 0(y ). 

Thus for large enough y the conclusions in the discussion above of the effects of the 

critical region on the /! root loci will not be valid, since they depend on the basic and 

critical region damping being of the same order in y. In particular, the behavior of 

the locus in which the root being destabilized by the critical region turns around and 

becomes more stable due to the eventual dominance of the basic damping is not possible 

2 

except for very small y, for which 0 (y ) effects are in fact negligible. Indeed, it was 
found in the numerical calculations that with y = 6 (y / 16 - 0. 375, i. e. , not verj' small), 
v near 1, and K p zero or small so the root is near 1/rev, that the jU locus does not 
turn around but rather eventually crosses into the RHP. The stability boundaries given 
in Fig. 9 are only valid then for truly small values of y/16. 


Order y Results 

O 

To carry the solution to order y“, it is first necessary to finish the order y 

solution. Considering only the case K - 0 and u £ n/2 for any n (i.e. , the basic 

R 0 

root), removal of the secular term from the equation for jS^ leaves 

^2 ill) (v + n) 




+ conjugate 


(69) 




v'T-i’ 





54£S£T‘ 
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The solution of this equation is 


s i =Re |W-- )e 


%*0 . V' M d 11 v o + “a - K p M e %K * n) 

' *«, ^ 2 e 

01 n^O n + 2 i/ q n 


( 70 ) 


and with 


i v ib 
0^0 




(71) 


K, 


*01 = *02 ( *2' • • > eX ’* 1 • S ’ 1 : ■'l + 2 M l + 1 M e 


(72) 


this completes the solution to y. 


The order y terms in the equation of motion give (dropping the common factor 


of y 2 ) 


2 2 2 2 

, ft , * ft , 2 ft , * - . 0« - 2 ft . 2 ft - 2ft>ft ft 

a^ o a0 1 


2^0 “ ^> 0 ^> 2 0 i o ' i^i 

- fr, 2 + *W V<V SW“i “ft ^ 


- |l- 2iv o^r + <_x i 2 - ■'i 2 - z v 2 + “/sV^I ° Xltl 


W. 


+ 1<- 2 Vi + M » - K p“e + 1 W '’n - 2 “-o -Sff 


_ i + Mo - K^M" .... 

- ft,, 2 ? °, p ■ °^° [-2\ i ft' 0 +n >- 2 Vi 


"02 *7* 2 

nA> n + 2 p q n 


+ M j 8 " K p M e + M /S 1 % + e 


Ml 


iP 0 V> 0 

e + conjugate (73) 


f 
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The secular term is, since v t n/2: 

i$ n . , e Vl L. *02 /. 2 2 0 lla 

ir-Vii=rnir| 2 ^^- + (\ + v i ^oWi)^ 


Regarding this as a differential equation for in terms of 0^, its secular term is 


with solution 


where 


^-V02 = ° 


^02 ^03 ( V * ‘ ) ® 


x i 2 + t -i 2 + 2 V 2 - M °s x i 
- 2it> 


2 0 2 
,M Al 


=iT 2 V2-^p“°e-(^M») -(J 


Thus the basic root, to 0(y ), is 


x = x 0 + yx i + y x 2 


‘V y [ l vH + i S^ M e 


xL _!i K M ° _(^if 

2V q [02 |/ 0 Ve \2V 0 M 0/ ' 2 ' J 


| M ? + it/ i+-^k p m°.-^(m| +(K p M 


2 0 

M e 
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2 2 

As reported above, there is no 0(y ) change in Re X, and the 0(y ) change in the fre- 
quency is dominant for small Kp, indeed for Kp = 0 the only change in the frequency 
2 

is 0(y ). To order M this root is 


X = - 


^ +i „ 

16 


K_ 

r 

2 

K 2 1 


1 + *- P - 

JL 

1 

l+(— 


16 2 

16 

. 2 

\ p / 


V 


2P 1 




which checks with the expansion to 0(y ) of the root from the small \x analysis (the 
hover root for 0(p.)). 

2 

The order y results would significantly alter the plots of the basic root loci 

2 

shown in Fig. 8. Extending the results for the critical region to 0(y ' would be much 
more involved because of the greater complexity of the solution for when 

v 0 = n/2. 

The Large y Case 


For the large y case, consider the general equation of motion, of the form: 
0 + v 0 = y [<M£ - K r Mq) ^ + (M^ - K p Mg) 0] 


(78) 


The small parameter in this case is the inverse of the Lock number. For y very large, 
the aerodynamics dominate the system. For y = ® the inertia and centrifugal spring 
terms (the LHS of Eq. 78) are negligible, leaving a first order system which does not 
depend on y, namely 


(M£ - K r Mq) 0 + (M^ - K p Mg) 0=0 


(79) 
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Reduction of the order of the equation of motion when the small parameter (1 /y) is set 
equal to zero is a characteristic of a boundary layer type of problem. The solution of 
the reduced equation is valid over most of the range in ip. As a first order equation 
however, its solution can involve only one free constant; thus it is not possible to start 
the solution from the two general initial conditions allowed for the original second order 
system. Furthermore, at certain points the solution of the reduced equation will exhibit 
singular behavior, indicating that the assumptions used to derive it must be reexamined. 
In general, there must be narrow regions in which the higher time derivatives are very 
large, so that inside the region the inertia terms are of the same order as the aero- 
dynamic terms and may not i>e neglected (other simplifications of the equation of motion 
are often possible though). If such a narrow region is used to connect a solution of the 
reduced equation to two initial conditions it is called a boundary layer ; if it is used to 
connect a solution of the reduced equation to another such solution on the other side of 
the layer, it is called a transition region . The solution of the reduced equation is called 
the main solution . More general terminology is inner and outer regions, and inner and 
outer solutions. Because the procedure for connecting the solutions in the inner and 
outer regions is central to the analysis of boundary layer problems, the entire analysis 
technique has been named the method of matched asymptotic expansions ; the name 
properly refers to the process of connecting the inner and outer solutions, but is usually 
used to include the entire analysis. Many techniques may be used to find the main solu- 
tions. Usually the technique used is a direct expansion of the dependent variable (j8) 
as a series in the small parameter (l/y). This technique is not satisfactory for the 
present problem because it does not yield a solution which is uniformly valid for all ip 





I 


r, 'T 



(which is required in order to find the roots); details of the application of this technique 
to problems are given in Ref. 3. It is also possible to use the method of multiple time 
scales to find the main solution. This technique is not entirely satisfactory either how- 
ever, since the equation obtained to lowest order is just the reduced equation given by 
setting l/y = 0 (Eq. 79). Since this is a first order equation, the solution gives only 
one root, which to lowest order (i. e. , for y = <») is independent of y. The reduction of 
the order of the equation means that one of the roots goes to -® as y goes to » ; i. e. , 
the solution corresponding to this root is exponentially small compared with the solution 
of the reduced equation. There is no way that the method of multiple time scales (as 

described here anyway) can find this root. The perturbation technique useful for finding 

tp 

both solutions in the outer regions is the use of a substitution of the form $ = exp J pd t/>, 
followed by an expansion of p(ip) as a series in the small parameter. 

The large y case is a boundary layer type of problem, which means that in general 
there will be several outer regions around the azimuth with a separate expression obtained 
for the solution in each region. Thus it will not in general be possible to find a single 
solution, uniformly valid for all </), from which the eigenvalues of the system may be 
found by inspection, as was possible for the case of small /i or small y. Instead it 
will be necessary to use the general techniques of the analysis of a system with periodic 
coefficients (as outlined in Appendix I). This entails obtaining the solution for one 
revolution of the rotor (one period), in pieces if necessary, each piece valid in a 
particular inner or outer region. First the main solutions must be obtained in the 
outer regions, but here it is necessary to find both main solutions rather than just the 
solution of the reduced equation. Next the method of matched asymptotic expansions is 
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used to connect the main solutions across the transition regions, or across boundary 
layers to initial or final conditions. Finally with the solution constructed over a com- 
plete period, the results of Floquet theory (see Appendix I) may be used to find the 
eigenvalues of the system from the initial and final values of the solution. 


Expansion in 


Consider first the equation with K = 0; write 

K. 


j9*exp J pd0 


p ’ yp -i tp 0 t ; p i + 


* 

P = p exp J pdi/) 

•• . 2 

P ~ (P + P ) exp J pdi/) 


Substituting these expressions for P, P , and j8 into the equation of motion, and collect- 

ip 

ing all terms of like order in y gives (exp J pd$ is a common factor in the entire 
equation, so drops out; the common factor of y n has also been dropped from the follow- 


ing equations): 


0(y 2 ): 


0(y): 2p_ x P 0 + P _! = P 0 + - K p M 0 


0(1): P Q 2 + 2p_ 1 P x + P 0 + V 2 * m £ P x 
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The order y equation gives p_ 1 = 0 or p_ l = M^, which give the two main 
solutions. 


First Solution 

2 

The order y equation gives p ^ = 0 so to order y find 



and to order 1 



Then the solution for )9 is 
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Then noting that 


exp 


; (ma)* 

- -Sn* 


M 


P 


= exp (- in Mp) = (M^) 


-1 


the solution for 0 is 




r p 




M a - KM, 


r r ph 

7 J M^d^ + J - d^ + 0(y ) 


P 


( 83 ) 


where P is a constant. 

A 


Eigenvalues 

From the continual appearance of in the denominator, it is evident that a 

transition region occurs where Mi = 0. This criterion means a transition region 

P 

occurs where the damping goes through zero. Alternatively, if is near zero, the 

P term in the reduced equation (Eq. 79) is much smaller than the P term, which implies 

that the inertia (0) terms must be included in order to obtain a differential equation with 

all terms of the same order; that is, there must be a transition region about the point 

where = 0. As it happens however, ($,M) is a negative quantity which never 

-1/3 

reaches zero; in fact, 8Mi < - (1 - 2 ) = -0.206 and even that value is never reached 

P 

unless p > 0.795; for p = 0, 8M1 =» -1 for all tb. Thus for the case considered 

P 

(K^ = 0), the main solutions are uniformly valid over the whole azimuth, and it is not 
necessary to deal with transition regions and boundary layers to find the complete 
solution. 
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With the solution over all ip, Floquet theory may now be used to find the eigen- 
values. Consider the first solution, to order 1. 


0 


-I 


M a - k p m « 


M 


3= JSj e 


0 


dip 


This may be written 


2ff 


M„ - KM, 


-*±/ 


0 = 0 ± e 


2ff "o M j3 


(84) 


where 


r Mo ” KpM« 2 ff M fl -KM 

f <*> - - / -Si / -S^- d * 

o p o p 


Using the fact that the aerodynamic coefficients M« Mi, and M. are periodic, it may 

p p V 

be established that f is also periodic in ip. Now Floquet theory states that the solution 
to a differential equation with periodic coefficients may be written in the form 


0- e^ u (0) 


(85) 


/ 7 \ 

\ 


f 

f 

*■ 


where 0 ^ is a constant, X is the eigenvalue, and u(t/>) is a periodic function (see 
Appendix I). Comparison of Eqs. 84 and 85 shows that the eigenvalues must be 


2ir 

2ir / 
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r m q - K P M e 
A ' Mi 


dip 









This result may be easily extended for both solutions to any order in y . Then the 


two main solutions give two roots: 


1 2 irJ MX * 

0 M 


+ ±j-l — e-. < -d^ + ofr*) 

y2ir I 


■*I 


f M -K M 


The symmetry of and means that 


2ir « 

f 

K 

Jn ► 


df/) a 0 


With this relation, the roots simplify to 
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and 


2ir 


2 it 



M^dj/) + K 


_L 
P 2n 


I 


M 


-M 


“d0 + 0(v _1 ) 


$ 


0 


0 


(87) 


{ 



/ 


Thus there are two real roots, one (X ) approaching -® as y increases to « 

(M^ is negative), the other (X^) approaching a constant; X^ is the root from the reduced 
equation. This behavior of the y root loci is expected from the small fi results; Fig. 1 
shows that for large enough y the locus is on the real axis, i.e. , there are two real 
roots, one approaching -® and the other -Kp for y -» ®. To lowest order X^ doet s ot 
depend on y, because it represents the balance of the aerodynamic damping and the 
aerodynamic spring only. The value of X^/(-Kp) for varying fi, and y = ®, is shown 
in Fig. 10; the movement shown takes place entirely on the real axis in the X plane. 

As for the small fi case (Fig. 1) the root is on the real axis, in the LHP if Kp > 0 
and in the RHP — unstable — if Kp < 0. The value of X^/ (-Kp) varies from 1 to 7/8 
for ft = 0 to «, with most of the change between fi = 0.5 and fi = 1; thus there is little 
variation of the root with fi (to 0(1)). The size of the 0(y *) term in X^ is indicated by 
the result for fi ■ 0, which is easily obtained (since the aerodynamic coefficients are 
constant then) as 


l l " " *P " y/16 


V 2 + K, 


+ o(y“ 2 ) 


This result agrees with an 0 (y *) expansion of the hover root frsra the small analysis. 

To lowest order X is: 

2 

1 2ff 

K 4 d * + 0(l).yM^."^2(-8M^) 

- w ■■ ■ ^ 
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The aerodynamic coefficient -8M^ is given in Fig. 7. For p £ 1 it has the value 
0 4 

-8M; - 1 + p /8; for large p it is asymptotic to 8/3 it p. This root becomes increas- 
P 

ingly negative as y increases, and also as p increases. The order 1 term in X 


is the negative of the lowest order (also 0(1)) term in X ; thus the behavior of this 

A 

term is also given by Fig. 10. 


Flap Rate Feedback 

When K = 0, there are no trims ition regions because Mi < 0 always. With flap 
R p 

rate feedback, K / 0, the same expressions for the main solutions are obtained except 
R 

that M£ is replaced by M£ - K^M^. The aerodynamic coefficient -(M^ • 
can become negative over regions of the disk for certain combinations of p and K^; 
i. e. , there may be negative damping over part of the azimuth range. When such regions 
of negative damping exist it means there must be transition regions about the points 
where the damping goes through zero. The main solutions obtained above are valid still 
in the outer regions, but in each region there are two constants, which must be matched 
through the inner region to the two constants of the next main solution. 



The criterion for the existence of transition regions is that there be negative damp- 
ing on some portion of the disk, i.e. , -(M^ - 

is usually positive, but may be negative on the retreating side for large enough p 

is too large positive, the negative values of Mg on the retreating 
side eventually dominate as p is inc reased, so there will be negative damping on 
the retreating side; if K R is too large negative, K^M^ eventually dominates on 
the advancing side and there will be negative damping there if p is large enough. 


(p > 0.641). If Kj^ 


KrMq) < Mjj * 8 always negative; Mg 
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Quantitative values of maximum and minimum K as a function of /J are given in 

K 

Fig. 11. For the cases with negative damping there will be transition regions (of 
“2/3 

width 0(y )) near where Mi - K M* = 0, which greatly complicates the analysis. 

p r y 

For these cases it is also expected that there will be other problems, including mate- 
rial computation problems, physical control problems, and large flapping amplitudes. 
The situation may be compared with stall flutter of a rotor in hover or forward flight, 
where a limit cycle oscillation is reached with the negative pitch damping in stall 
balanced by the positive damping below stall, resulting in high amplitude pitch motions 
and large control loads. Thus while a region of negative damping does not necessarily 
mean there is a flapping instability, it does mean that there are many problems — 
analytical, computational, and physical — so requiring -(M^ - K R M^) > 0 is a reason- 
able design criterion. This criterion provides a maximum and minimum K for a 

R 

given /l. The limits of K from this rule are much easier to obtain than actual stability 

R 

boundaries; and Fig. 11 shows that although conservative, it is not a serious restriction 

for /I less than 1 or 2. For large /i it is a serious limitation (for large /i, 

K r = 2/3/1 (1 + 7/12/1) and * - 2/3/1 (1 - 7/12/1)), indicating that M Q 

max min 

(blade pitch) is not very good for flapping rate feedback then. Time varying K might 
work better but it would have to be programmed with /I probably. Although the deriva- 
tion of this rule has been based on the large y case, the criterion of no negative damp- 
ing has nothing to do with y, and so should be a reasonable criterion for all y. Indeed, 

the criterion K > -1 for /i - 0 is the same as from the small /i case, where it is a 
R 

true stability criterion, and valid for all y. 
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The Large M Case 





For the large y case the general flapping equation of motion, as given in Eq. 1, 
is used; the small parameter in this case is the inverse of the advance ratio. For jU 

very large, the aerodynamics again dominate the system. When y goes to infinity, 

• %* 

the j3 and & terms are arbitrarily large compared with the (9 term because of the 

influence of y in the aerodynamic coefficients. Thus this problem is also of boundary 
layer type, and the solution is sought as for the large y case in terms of outer solutions 
and transition re^'cns. When ^ 0, the aerodynamic damping, K^M^, is the same 
order in y as the aerodynamic spring term, M- - K_M-, (namely 0(/i 2 ), see Eq. 1), 

p P U 

so setting l/y to zero reduces the order of the system. The reduced equation gives 

one main solution, and the other will be exponentially smaller (or larger). When 

K = 0 however, the aerodynamic damping. Ml, is of 0(y) while the aerodynamic 
R P 

2 

spring, M fl - KM is 0 (y ) (Eq. 1), so in order to obtain an equation for the outer 
p P o 

solution with the proper ordering of terms it is necessary to include the inertia term 
(/9) even in the equation for the outer region. That is, for y -* <» the aerodynamic 
spring must be balanced by the inertial forces, which leads to an equation of the form 


j8+fiC#) 0 + ATKM j8 = 0 


( 88 ) 


The solution of this equation is either a rapid sinusoidal oscillation with frequency of 
0 <M), or a sum of exponentials with time constants of 0(/i _1 ), depending on whether 
the aerodynamic spring is negative or positive (the criterion is a bit more complicated 
really, but that statement will do for the present discussion). Now the aerodynamic 
spring changes sign in the middle of the advancing side and again in the middle of the 
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- 2/3 

retreating side, and at each point there is a transition region (of width 0 (/x '")) across 

which the solutions must be matched. There are also transition regions (of width 
-2/3 

OQl )) between the advancing and retreating sides of the disk (i. e. , around ip = 0 and 

180°), through which the main solutions must be matched. Thus in contrast with the 

large y case, for large JL4 it is not possible to find an outer solution uniformly valid 

for all tp. Rather it is necessary to go through the entire procedure of matching the 

main solutions through the transition regions and by that process construct a solution 

over one rotor revolution. Then the results of Floquet theory (Appendix I) may be used 

to obtain the eigenvalues from the initial and final values of the solution. The procedure 

for finding the main solutions will again be based on a substitution of the form 
jp 

P = exp J pd ip with p now expanded as a series in /i . The matching techniques of 
the method of matched asymptotic expansions will be illustrated in the treatment of the 
transition regions. The procedures required here are reasonably straight forward, but 
in general the matching techniques can be quite complicated, particularly when solutions 
are sought to higher order. The reader is directed to Ref. 3 for more details of the 
method of matched asymptotic expansions. 


In regions (1) and (iii) of the rotor disk the differential equation has the form 


•• o 

j3 + v j3 = - r y 


f + sin»j!) + ^ + j M sin 0 + ~ (H sin 

M cos ^6 + 4 ** Bin ty + K p (g + 3 M 9ta $ + 4 0* sin ^) 2 ) 


(3 


(89) 


where r is a constant with the value +1 in region (i) and -1 in region (iii). Recall 
from the discussion of Eq. 1 that region (i) is the advancing side of the disk, where the 




T 


/ . 


blade has normal flow over its entire span, and region (iii) is the range of if) on the 


retreating side where the blade has reverse flow over its entire span. As n - co, 

region (iii) occupies nearly all the retreating side, with the exception of 0 (jz *) bands 

near if) = 180° and if) = 360°. In Eq. 89, M appears in the aerodynamic coefficients 

nearly always in the combination M sini/); any assumptions made about the order of terms 

based on the order of fx will be violated then if sin 0 is small enough; this is the origin 

of the transition regions near if) = 0 and 180°. When K = 0, there are also transition 

R 

regions in the middle of the advancing and retreating sides around the point where the 

aerodynamic spring goes through zero. These regions arise in the analysis because the 

aerodynamic spring being zero or very small will again violate the assumptions made 

about the order of the terms; physically they arise because there must be a transition 

between the solutions on either side since they have very different behavior, namely 

sinusoidal oscillation and exponential decay or growth. If K ^ 0, such transition 

K 

regions are not required; only the transition regions between tb lvancing and retreat- 
ing sides are needed. The case with K = 0 will be examined firbt. 

R 

Expansion in U 

Consider the outer regions, where fi sin if) is of order fi and all time derivatives 
of 0 are of the same order; this means the regions (i) and (iii), away from the boundaries 
near if) = 0 and 180°. The equation of motion in the outer region is then Eq. 89 with 
K = 0. For K * 0, the analysis is simplified if first the 0 III) aerodynamic damping 
is removed from the equation of motion (which is of the form of Eq. 88 to the lowest 
order); this is accomplished by the following substitution 


r fi cos if) 

0 = e y (90) 
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With this substitution, Eq. 89 for K = 0 becomes 

R 


• • 2 y 

y + V y = - r J 


y-ryy “ M ^ |sin(|)|^ + sin^ 

+ ju cos0^ + j M sini/J) + K p ^ + sin0 + j (/i sinipf^j 


( 91 ) 


Then the main solutions are found by use of the substitution 


* 


y = expj pd0 


(92) 


with p expanded in a series in /jl 


-1 


P = MP.i + P 0 +J p i + -.. 


Main Solutions 

Substituting for y and collecting terms of like order, the equation of motion gives 
0(M 2 ) : ^ | sin ip | ^ | sin ip | + cos ip + K p sin ip) 

0(M): P_ x + 2p_j P 0 = “ r g P _ x “ r y (- ^1 sin# | + ^ cos ip + j K p sin^ 

2 

From the order fl equation, one obtains: 

P^ - ±*/^|sin 7 1^1 sin 0 1 - cos </) - K p sin i/)) (93) 

jp 

The solution for 0 has a factor of the form exp fJ, J p_ 1 d0; the double sign in p gives 
the two main solutions. When the quantity under the square root is positive, p is real 
and there results main solutions with exponential decay or growth (with time constants 
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of 0(|U *)); and when the quantity under the square root is negative, p ^ is imaginary 
and there results main solutions with sinusoidal oscillatory behavior (with frequency 
of 0(/i)). When the quantity under the square root is zero, there is a transition region. 
On the advancing side (region (i)), p ^ is zero at 



(94) 


and on the retreating side (region (iii)), 


0 = tan 

o 



(95) 


So there is a high frequency oscillation on the rear of the disk, exponential solutions on 
the front, and transition regions between the two types of behavior, p is also zero 
when sin0 = 0, i.e. , at 0 = 0 or 180°, so transition regions will also be required near 
the edges of regions (i) and (iii). In addition there is region (ii), which for large H is 
an order (1 * small band on the retreating side near 0 = 0 and 180°; in connecting the 
main solution from the advancing side to the retreating side it is necessary to go through 
this region as well as through the transition regions. In general, whenever p is 
small the assumptions made about the order of the terms in obtaining Eq. 92 are violated, 
so the main solution can no longer be valid there. This criterion gives the four transi- 
tion regions. 


The order /i equation gives 


eo" 


2p 


-1 
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from which, using the relation 


>'p.i 




one obtains 

Jp 


J 


V* - 


-2*1 P -l! r 2 


r 


’ 9^l Sin ^l + 12 C0S ^ + 3 K p Sin ^ 

_ — “ “ 


dip 


With the above expressions for p ^ and p^, the solution for ft (to 0(1) in p) is 

obtained by using the substitutions for ft and y, i. e. , Eqs. 90 and 92. There are main 

solutions in four regions, which will be called quadrants (although they are not really 

so, since ipf. and ip^ are not equal to 90° and 270°). Each quadrant is bounded by 
1 3 

transition regions ; the ranges of the four quadrants and the main solution valid in each 
are as follows. 


1st quadrant: 0 < ip< ip 


t. 


i * 






+ conjugate 
2nd quadrant: ip t < ip< n 


(C, 


/ Jp J> g \ 

MM J /Tdip-J 'fed ip) 

+ i C„) e ' * *1 ' 


(96) 




0 0 _ 

M J ST dip - J* yfdip 

C.e ** *2 


~ r ‘^5sX-- 


+ C. e 
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0 0 g 

- M J /Td^ + J -»d(l) 

h 0 2 


(97) 
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3rd quadrant: it < ip < ip t 


- u cos ib + ib , 
j8 = e 12 16 (ff l/4 


0 0 g 

Mj vfd^ + J jfdip 

C e ^ *B 


ip 


+ C e 
6 


-Mj ./fdj/i-J^ ^jdtf 

03 


r 

» , 

03 


(98) 


4th quadrant: ip t < ip < 2ir 

l 3 


- u ^ — cos ib + ib / 

S = e 12 16 Hf 1/4 


+ conjugate 


(C ? + i C g ) e 


>P >P cr 

i( M J /=*d* + J ^=dd 
V 0 4 >P/ ' 


(99) 


where 


f = J | sin ip |^r | sin 0 1 - cos ip - K p sin ipj 
g " 2(12 008 0 + 3 K P Sin 0 " 9? |8ln *•) 


and Cj, Cg, C g are constants. The matching procedures will result in connection 

formulas through the transition regions, which will give the two constants of one main 
solution in terms of the two constants of the main solution on the other side of the 
transition region. In addition, the constants will be matched to arbitrary initial condi- 
tions at a certain point on the disk. The quantities ip , ip 0 , ip , and ip. in Eqs. 96 to 99 
are also constants, which must be in the appropriate quadrant; they are not free 
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constants since a change in them must be accompanied by a change in the value of the 
C's. These angles may be given any value convenient to the analysis; it is most con- 
venient here to leave them arbitrary since they will drop out of the final result anyway. 


Transition Regions 

Consider the transition region near if) (meaning in region (i) or ^ in 

region (iii)). This transition region has a width of The reader is directed 

to Refs. 3 and 4 for illustrations of methods for finding the proper width of a transition 

region or boundary layer. The technique involves assuming d 0 = 0(M °); then $ is of 
_2n 

order fi , and similarly the order of all terms in the equation of motion may be 

found in terms of n. The exponent n is determined from the criteria that the resulting 

equation, to lowest order in fi , must a) include the highest time derivative (j3) and 

b) must produce solutions capable of being matched to the outer solutions. It also helps 

-2/3 

to know what to expect of certain types of problems; for example, a width of 0(p ) 

is typical of transition regions for equations of the form of Eq. 88. 

-2/3 4/3 

Assuming d$ = 0(p ), Eq. 91 becones to lowest order (0 (p. )): 


•• , ■ / 2 5f 

y-y (♦-*,) 4 sj 


( 100 ) 


where 





2 1/3 

Alternatively, write z = (r (y/i) fi ) (ip - tpj, substitute for ip in the differential 
equation for y, and then obtain to lowest order the equation 



z y = 0 


( 101 ) 


This equation is a standard form, the solutions of which are called Airy function. There 
are two independent solutions, denoted Ai(z) and Bi(z). These functions may be written 
in terms of Bessel functions; however, the general behavior of the solution in the transi- 
tion region is not of interest here. Rather the solution in the transition region is only 
to be used to find a connection formula between the neighboring main solutions. For 
this purpose all that is required is the behavior of the solution for very large z. 

Writing the solution of Eq. 101 as 

y = z/rTa Ai(z) + v^T b Bi(z) 


where a and b are constants, then the behavior for large z is: 


-1/4 < . C v 

z-*o>: y~z (ae +be) 


z -* y ~ |z| 


r n 

-1/4 /b \ J 4 K 

- i a) e e + conjugate 


where 


c-| w s/ 2 -f«.*-v 8/2 


The matching procedure consists of finding the limit of the outer solution as 
ip -* ip , and the limit of the inner solution as z -* A<», and requi ring that the two limits 
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have identical behavior. This criterion gives the constants in the inner region in terms 

of the constants in the outer region. The matching procedure is considerably more 

complex if higher order solutions are involved. Consider first matching from the first 

quadrant to the second quadrant, through the transition region at ip = lp t . The outer 

l l 

solution in the first quadrant is given in Eq. 96. As 0 -* 

$ 4> t ip 

J* /Td</> = J /Tdip + J SJdtp 


-I'r- 


so the main solution for ip ^ is 

, , 1/4 

y outer' - !l * ~ VI / (C l + lC 2 )e ' 1 *1 /e 

+ conjugate 

To the inner solution, the outer region in the first quadrant appears as the limit 


z -* -oo, and in this limit 


-(y)' lA 2 l*-* t f 1 / 4 (f +i ») 


- 1 ” 


e e 


+ conjugate 


Then requiring (y . ) . , 3 (y. ) 

w outer inner 


gives 
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and the inner solution becomes, for z -* ® 


V 

'inner 



( - Mfl *-*/ 2 

la e 


+ b e 


/2 

3 



and the matching criterion gives 


Mj *vTd*-J*7.d» 

e *2 V 



C 


4 



st 

/Td^ + J yj# 

^2 



(103) 


Next combining the results of matching between the first quadrant and the inner solution 
{Eq. 102), and between the inner solution and the second quadrant (Eq. 103), gives the 
connection formula between the first and second quadrants: 



/ . 



(C + i c 0 ) e 
' 1 2' 


\ *1 *1 /. e 4 


= e 4 Ac e *2 

2 3 


<l> t 4> t n 

nj V!d*-; '4d* 


<l> t be 

-(if v/fd^ + J 7?d0 

+ «C e *2 *2 

4 


Similarly, the matching procedure on the retreating side around ip t gives the connection 

l 3 


formula between the third and fourth quadrants. 

/ ^ h s \ r 
l W A - 7d < ,+ J\ 73 d * if , 

(C 7+ lC 8 )e\*4 *4 7 =e 4 f < 




TCe 
2 6 


•iC e *3 
5 


4. A g - 

/Td^i + J ysdil> 

, h h 


Now consider the transition regions near |(iaOor 180°. The main solutions indi- 
cate that there must be transition regions at the edges of regions (i) and (iii), but there 
is als'' "ogion (ii) in between regions (i) and (iii). The extent of the three regions is 
defined by (see Eq. 1): 

region (i) <• sin > 0 

region (ii) - 1 < M sin ^ < 0 
region (iii) fi sin $ < - 1 
-2/3 

The transition region has a width of order M ; thus region (U), which has a width 
of only order fi lies entirely within the transition region. Region (ii) appears then 
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J* 


as an interior region which has no effect on the solution to lowest order. This means it 
is not necessary to use the more complicated aerodynamic coefficients of region (ii). 

4 

It also means that the appearance of p terms in the aerodynamic coefficients is 

deceptive, since these do not appear except in region (ii), in which region (p sin 0) n 

is of 0(1) or smaller for all n. The true order in p is given by the aerodynamic 

2 

coefficients in regions (i) aid (iii), and these are of 0(p ) at most, as expected of 
aerodynamic forces. 


The proper matching procedure is to find the solutions in the transition regions 

at the edges of regions (i) and (iii); in addition, the solution is found in the interior region 

including region (ii) and the neighboring parts of regions (i) and (iii) where p sin 0 

is of 0(1). Then the matching process proceeds from the main solution in region (i), 

to the transition region at the edge of region (i), tc the interior region, to the second 

transition region in the edge of region (ii), and finally to the main solution in region (iii); 

by this process the connection formula between regions (i) and (iii) is established. With 

the substitution x ■ p sin 0 (x is assumed to be of 0(1) so 00 - 0(p 1 )) the equation of 

2 

motion for the interior region is found to be, to lowest order (0 (p )): 


the solution of which is 



*- 0 




(106) 


where ^ and 0 2 are constants. Thus matching through the interior region is just a 
matter of matching the displacement and slope of the neighboring transition regions. 
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-2/3 -1/3 

To the transition region (of width 0(/L l )) the interior region appears 0(p ) small; 

-l/3 

then in terms of the transition region variable z = fl x, tne matching process is to 
be carried out at the limit z -» 0. 

2/3 -2/3 

With the substitution z = /x sin $ (which implies an 0(p. ) region around 

if) = 0 or 180°), the equation of motion for regions (i) and (iii) (Eq. 89) becomes, to 

4/3 

lowest order in ju (0(p )): 


^-f±£|z|0 = O (107) 

dz 


The plus sign applies on the back of the disk and the minus sign on the front. In terms 

of the variable z, the ad mcing side is given by z > 0; and the retreating side by 

-l/3 

z < 0; on f h-i sr - region (ii) appears as a negligibly small (0(p )) area at z - 0. 


Consider the back of the disk, i. e. , from the fourth to the first quadrant. The 
differential equation for the transition region is (Eq. 107) 

d -^|zlA = o 

dz 


The solution again involves Airy functions, and may be written 
|» At[-(J) 1/3 zj + b Bi (*) 1/3 zj 
|a* Al[(*f /3 zj + b* B)[(*) 1/3 zj 


where a, b, a*, and b* are constants. Matching the displacement and slope at 
z = 0 gives a* *\/3 b and b* * aA/3, so 
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yS =-s/3 b Ai [(^) 1//3 z 



z < 0 


The asymptotic behavior of the solution in the transition region is chen 

- 1/2 


z -*■ oo: 


ft- 


z -* - oo: & ' 


iu/ z J 1 


-[© 1/3 « j ' 1/4 -'' 1/2 


' \ 
b - i a i 4 i£ 

— “ — I e e + conjugate 


/ \ .£ 

a-i3b\ 1 4iC, 

CUT") 6 6 + conjugate 


3/2 

where £ = \/y/3 z , note that the solution has oscillatory behavior, which is the 
proper behavior for matching the main solutions on the rear of the disk. Now the main 
solution in the fourth quadrant (Eq. 99) becomes, as # -* 2 n 


{$ -* e 

outer 


W 2 ’)(I m -2/3 |s ,)-1/4 


+ i -J_) e V "4 *4^ / 3 + 


< C 7 + 1 J 8> 


conjugate 


and the main solution in the first quadrant (Eq. 96) becomes, as i/i - 0 


/3 -* e 

outer 


" 10 1 C 


12 T 16 2 ”/y -2/3 , \-l/4 

(4" ") 


4 f d*\ + l^|2| 3/2 
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conjugate 
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Then the matching process, requiring ^ ou t er )0_2 ff = ^inner^z- <» in ^ ^ ourtl1 y ua drant 

and (j3 ) = (/3 ) , „ in the first quadrant, gives the connection formula: 

inner z=® outer!/)=0 


g 


(C^ + i C 2 ) e 


imT ^ d *-J 7^ m z_^l 


h 


4n 


= 2(C 7 + iC 8 )e 


*1 / e « 16 

27T 2ff 

i(/lj /Tdi^ + J 7%dl/)l 

1 h v 




2tt 


-i mJ* /-Td^ + J 


i (C ? - i C g ) e 


(108) 


Finally consider the front of the disk, i. e. , from the second to the third quadrant. 
The differential equa‘ i >n is (Eq. 107) : 


-4 |z| 3 = 0 

dz 


(109) 


The solution is 


3 = 




z > 0 


z < 0 


Matching the displacement and slope at z = 0 gives a* =/3 b and b v = a//3, so 


^ -/3 b A1 (-^) 1/3 *) *% B1 z) 


z < 0 






1 ^T\ - 


"T 
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The asymptotic behavior of this solution is 


z - >: 0 1/4 77-1/2 [f e "^ + b e ^j 

ft //1V/ 3 , W 4 ~- 1/2 RIk . _a_ cl 

*~\u] ,z y |y + 7 f e J 


where C =/y/3 izi ; this solution has the exponential behavior required for matching 
on the front of the disk. The main solution in the second quadrant (Eq. 97) becomes, 


as 4> - 


/3 -* e 

outer 




M / /Tdlfi-J - Jj dip Kl 3 ^ 2 


-(*/ yid**.” ^d* + ^| S , 3/2 


4 C. e 
4 


Matching this to the inner solution, i.e. , requiring (/5 . ), ■ (A ) , gives 

outer <p=7i inner z™ 


Mj /Td^-J 4d^> 




JL/W/6 a 

yft \ 4 u/ 2 


-mJ* /*d* + J* 4d0 _ U ^L_X 


M 12 ■ 16 3ff 


- J-fc-Y 
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The main solution in the third quadrant (Eq. 98) becomes, as I 


a ^12 + 16 n (X n -2/3 \-l/4 

P . -e (TM izi 

outer \4 / 


* _ , / g . yv 3/2 

/? dl/) + J yjdi/) + -^-izi 


/ _ r ff g . ^ 3/2 

- Mj /fdl^-J yydjJ)-— tzi 

+ C e *3 *3 ‘ 

6 


Matching this to the inner solution, i. e. , requiring (^ inner ) z= _ co ~ ^ 0 uter^=ff’ giveS 


"I, /Td* + J Tfd# JL + JL, 


C e *3 


12 16 _ A./X.V/ 6 JL 


-Mj /Td 0 -J 4 d«/) u 2 _ + 2 _., 

C e ^3 e 12 16 " = ^V /6 ^ b 

C 6 e /i\4M/ 2 ° 


Then combining Eq. 110 and Eq. Ill gives the connection formulae: 


Mj 4d* _ U I_^- 


6 " 16 " /I 


= fC e 
2 5 


6 

Mj /Fd^ + J* "/r dif) 

h h 


-Mj /Fd«/) + J d«/> 


6 16 


2 „ 

/3 "6® 


-Mj* /Fd^-J y5d0 
*3 *3 
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The derivation of the connection formulae across the four transition regions completes 

the construction of the solution around the disk. It is also necessary however to start 

and finish the soiution with initial and final values of jS and 0 at some point. It is most 

convenient to start and finish the solution at 0 = ff; this is the middle of the interior 

region, where it is easier to match the two constants in the solution to arbitrary values 

of 0 and 0 than at any other point on the disk (to lowest order, the solution to Eq. i06 

in the interior region is linear in x, so the two constants are easily related to 0 and 
• 

0). In terms of the transition region, the slope and magnitude of the solution are to be 
matched to 0 and 0 at z = 0. Now the solution of Eq. 109 becomes, for small z 


from which 


0 = d x (a+/3b) d 2 (- a+/3 b) z 


0(ff) =d x (a+/3b) 

b m d 2 (a -/Jb) 


d, d n \ 4 , 


,2\-l/3 


2/3 b = 


d l d 2 ' 4 ' 


,2\-l/3 


where dj and d g are constants associated with the Airy functions (dj - 0.355, 
d ? - 0.259). Then for initial conditions, Eq. Ill relating C g and Cg to a and b 
gives the main solution in the third quadrant in terms of initial values of 0(ff) and 0(ir). 
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For ending the solution, Eq. 110 relating C and C to a and b gives the final values 

o 4 

* 

of £(ir) and $ (ir) in terms of the main solution in the second quadrant. 


Eigenvalues 

With the solution in each of the four outer regions, and the connection formulas 
between them, the complete solution may be constructed over one rotor revolution. 
Floquet theory for a single degree of freedom, second order equation (see Appendix I) 
shows that the eigenvalues are given by the quadratic equation 

(e^V - «3 r + Bp) (e* 2 ") + B R B p - B R B p =■ 0 (114) 


where 


0 p = j9(ir + T) 
j3 p = /3(jr + T) 
fl R = 0(ff + T) 
|3 R = /3(ir + T) 


due to ]3(rr) = l, 
due to 0 (ir) * 1, 
due to 0 (ir) * 0, 
due to 0 (ir) = 0, 


0(ff) = O 


kit) = 0 


kit) - 1 

kn) = 1 


and T is the period of the system (T = 2 n here). 


Combining the connection formulae (Eq. 104, 105, 108, and 112) and the initial 
and final value formulae (Eqs. 110, 111. and 113) results in the following expression 
for 0^ and ^ at $ = n + T in terms of 0^ and 0^ at $ * it: 




_ — 


T— 
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(115) 


where 


F = M J^ 1 /3 cU/) - # 

1 0 0 


IT 


IT 


g 


F 9 = mJ yjdj/) 


<lv 




A 


Fg = M /^ 3 /I di/) + f d0 


IT 


2 n 2ir g 

F 4 =mJ /^d0 + J yq? d„ 


*L 
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Obtaining /3p, j9 p , jS^, and ^ from Eq. 115 and substituting into Eq. 114, there 
results (after some manipulation) the following equation for the eigenvalues: 




(117) 


where 


. Ijl -< F 2 ' V 

b *2 4® 






[2 cos (F x + F 4 ) - cos (F x - F 4 )J 

4 ^2 ^ ^*3^ r 1 ) 

|e 3 [2 cos (F x + F 4 ) + cos ^ - F 4 )J| 


( 118 ) 


l- 


-«A. 

, /? 
. Vt 

A 


% 

K- 

A 

m 




•-«> 



So b is a function of /x, y, and K p . The solution for the roots is then 



y 1 -1 

- u ± — cosh b + ni for b > 1 

677 277 

y 1 -1 

- a 7 — ± i — cos b + ni for -1 < b < 1 

677 277 

- u ± ~ cosh * Ibl + ~ + ni for b < -1 

” 677 277 1 I 2 


(119) 


where n is some integer. This result shows the typical behavior of the roots of 
periodic systems. For b < 1 the damping is fixed at - M (y/ 677 ) with a change due to 
b in the frequency; for b > 1 the frequency is fixed at n/rev with a positive and nega- 
tive change due to b in the damping; for b < -1 the frequency is fixed at n 4 /rev with 
a positive and negative change in the damping. The critical region boundaries are 
given by b = 1 and b = - 1 . 


The general character of the critical regions and instability boundaries in the 
y-M plane, as obtained from the solution of Eq. 117, is sketched in Fig. 12. Because 
M is large, it happens that |b| is much greater than 1 almost always, so the critical 
regions dominate the behavior of the roots. Because of the cosine terms in b, the 
sign of b changes regularly; b must ox course go through zero then, but it does so 
very quickly, so there is only a very narrow band between the ImX- n/rev and the 
Im X = n + 4/rev regions in which |b| < 1 . When |b| < 1 , the real part of X is -M (y/oir), 
i. e. , the root is stable for all M and y ; thus there must always be a band of stability 
surrounding the transition from n/rev to n 4/rev. These characteristics are illus- 
trated in Fig. 12. The locus between tlv critical regions has a rather fine structure 
which would be difficult to obtain numerically. A root locus for varying n or y 


■ — r ^ V 



(a vertical or horizontal section in Fig. 12) in the vicinity of a critical region boundary 
would in quick succession move from the RHP (unstable) to the LHP (stable) with fre- 
quency fixed at n/rev, rapidly move from Im X = n/rev to Im X = n + £/rev in the 
RHP with damping given by - M (y/6ir) (which would be nearly constant because the 
critical region boundaries are so close) and then moves from the LHP into the RHP 
with frequency fixed at n + £/rev. 

Figure 12 shows that for a given fi the system is stable for a large enough y. 
Positive K p is stabilizing, tending to decrease the size of the instability regions; 
negative is destabilizing in this sense. The rotating natural frequency of the flap 
motion, u, does not enter the high M case to order p^ (the aerodynamic spring 
dominates the centrifugal spring until order p^); this is consistent with the fact that 
the critical regions dominate the high behavior, so the frequency of the motion is fixed 
at a multiple of i/rev. 

A comparison of these analytical results with the results of numerical calculations 
indicates that the high JU solution is good down to /J = 2. 5 or so. Thus numerical calcu- 
lations are required to join the loci from fi = 0. 5 to 2. 5 say (for y neither small nor 
large). The behavior theoretically predicted for the locus at large fi (in particular the 
rapid movements between Im X * n/rev and n+^/rev, and perhaps — for y not too 
large — between the RHP and the LHP) actually does show up in the numerical calcula- 
f ions of the stability (above M 3 3.0 say); such behavior of a numerical solution might 
be questioned without the perturbation solution to provide a guide to what to expect. It 
is unfortunate that the boundary of the instability region for y/16 of order 1 is first 
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encountered at moderate jU (around JU = 2.25 for small K p ; see Fig. 12) and so cannot 

be obtained by perturbation techniques (to the order explored anyway). Because of the 

small time constant in *he main solutions (0 (jl 1 )) and the four transition regions (of 
-2/3 

width Q(ji )), a numerical calculation of the roots for truly large would be difficult; 
the perturbation theory handles these singular problems analytically, and the calcula- 
tions that remain are nonsingular, short, and simple. 


Flap Rate Feedback 

The use c. , llap rate feedback (K^ / 0) changes the solution for large fi funda- 
mentally, because the aerodynamic damping (-K^M^) is then the same order as tbr 

aerodynamic spring (M« 

P 

then, and only the transition regions near if) = 0 and 180° are required. The main 
solutions are obtained using the substitution 

J) 

3 = expJ p dif) 

with p expanded as a series in /i: 

p. ( . 2 p_ 2 + 4p_, + p 0 +ipi + ... 

Making this substitution in Eq. 89, the terms of like order in give 
00* 4 >- p 2 2 .-ry^K H i(»to*) 2 P _ 2 ] 

°(K 3 >: 2p.jP.! ■ - r y .1. * (£ * r. 2 + K R ± (sin « 2 p.,] 

0(P 2 ) : p 2 ! +2p.j + P. 2 - - r y[i 0 + K r ) P. 2 +; sin *(j+K R ) p -l 

+ Kj! j (sin W 2 P 0 + j sin #(ooa * + Kj, sin #)1 

f -T- — — — 


KpMg)* derivation of the main solutions is simpler 
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OW): 2p_ t P 0 + 2p_ 2 Pj + p_ x - - r y [i d + Rjj) P j + i Sin *(| + Kjj) p 0 

+ Kr 4 ( sin ty) 2 Pj + c °s 0 + J K p sin 

The first solution is given by p_ 2 * 0. Then the order n 3 equation gives also 
2 

p_j =0, and the order \i equation gives 


P 0 = 'K~ ( cot ^ + V 
R 


or 


J ! S 

J P Q d0 = - ~ Xn sin 0 - ~ 0 


I 

t 

v 


The order n equation gives 


% 


P, - 


2 L . 1 \ cos £ + 2K P l 


1 ' 3K H'( 1+K J(.in*) 2 ''3K n i!81 "* 

XV 


or 


J 1 3K r \ A Kjsini/) 


2K p * 

+ 2 intan f 


3K r 


T hen the solution for p is 


V V " 


3Kp 


X. n l6 M2K R M»M 3K V K ; 

/9 » jS, (sin 0) e " (tan*) 11 e R R 


K R “r 


3 K r l Kj 


+ 0(M -2 ) 


(12 


- 





\ 

/ 


where 0^ is a constant. This is the solution of the reduced equation, and so is independ- 
ent of y to this order, i.e. , is the result of a balance of aerodynamic damping and 
aerodynamic spring terms only. The solution is not valid when sin 0 is too small; 

and it is exponentially growing (unstable, at le> .st in the outer region) if K /K < 0. 

P R 


The second solution is given by p = - r (y/4) K (sin 0) or 

*2 R 


/ P_ 2 = - r 4 K r ^ * " 4 8ir * 2 


The order fj, equation gives 


r g (1 + 2 Kjj) sin 0 


r* 1 


J P_J W ■ r Q {I + 2K R ) 008 0 


and the order fi equation gives 


1 - 2K 


p o- r 8< 1+K R> + nc ico ** + ir 


> i - 2K » rs v i 

I Pq ^ ^ + |JSl ” r 8 ^ + K Rj ^ 


Then the solution fcr $ is 


0 * 0 2 (sin 0 ) R exp j- r jj[m 2 Kj. (20 ~ »ia 20) - M J (1 + 21^) cos 0 

♦ 2(l + K R )0j + j^0 + O(M’ 1 )| 

— 1 ; —t- — — — ’ 





I 


) 
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where 0 is a constant. This solution does involve y, i.e. , it involves the inertia 

Lt 

terms in the equation of motion; this solution is also not valid when sin 0 is too smell. 

There are transition regions between the advancing and retreating sides, near 

0 = 0 and 180°. As for the K =0 case, region (ii) is an interior region lying entirely 
-2/3 

within the 0 (U ) transition region; the solution is matched through th/3 interior 

region again by matching ine displacement and slope of the transition region solutions 

2/3 

at z = 0. With the substitution z = M sin 0, the equation of motion in the transition 

4/3 

region becomes, to lowest order in (Q(H )) 



*df +(i ) =0 


( 122 , 


, where the plus sign applies on the back of the disk and the minus sign on the front. 

| It can be shown that the solutions of this r quation will have tfcj proper asymptotic- 

behavior for matching to the two main solutions. Unfortunately the solution of this 
equation is not available in terms of classical functions. The behavior of the solutions 
for large and small z would have to be found, probably largely by numeiical methods, 
before the matching procedure could be ca. ied out. 


i 






-r 






Applicability of the Four Cases 

This section will consider the ranges of fi and y over which the four 
cases investigated above are useful. Perturbation theory is based on the 
expansion of quantities in terms of a very small or very large parameter. 
In r?any problems however, the results are useful, even quite accurate, 
far beyond the limits for which they are theoretically valid. It is these 
perturbation solutions, extendable up or down from truly small or large 
values of the perturbation parameter, which are of most value. They 
may be found only by comparison with exact solutions, usually obtained 
by numerical methods, for moderate values of the perturbation parameter. 
Another question about the range of validity of the solutions arises in this 
problem because there are two parameters, fi and 7 , which are avail- 
able as perturbation parameters. In the perturbation analysis based on 
one parameter, say y small or large, the solution is derived under the 
assumption that the other parameter is of order 1, e. g., fi = 0(1). This 
raises the question of the validity of the solution when the other parameter 
is itself very small or veiy large. It may happen that the results are 
still valid when the other parameter is outside its assumed range, but 
this must be checked in each case. The ordering process of perturbation 
techniques provides a quantitative framework for making this check. It 
is usually quite simple to determine what range the other parameter must 
have so that the assumptions made about the order of terms are still valid. 


O*. 


T 
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Extending the analysis into ranges of the other parameter where the ordering 
assumptions were violated is another matter; it of course means an entirely 
new case to be considered and analyzed by perturbation techniques. The 
results of the analytic solutions obtained above were compared with numerical 
calculations (performed by the author) of the roots of the flapping equation, 
primarily for moderate and small 7 , over a wide range of fi ; these calcula- 
I tions were similar to those reported in Ref. 2. 

! 

The small (i results are good out to H =0.5, which is a very useful 
range. These results are valid for all y , since the order of y does not 
change the terms retained in the small H analysis. The large ft results are 
valid above = 2. 5 or so, which is also a good range. Here however y 
either very small or very large will violate the assumptions made about 
the order of the terms in deriving the large ji solutions and so the results 
may not be valid in these corners of the y-fi plane. 

The small and large 7 results, to the order investigated, are really 
useful only for truly small or truly large y , although the results are quite 

informative. The small y results are accurate up to y = 2 or 3. The 

i 

| limitations of the 0(?) solution, as discussed in the small y analysis, pre- 

vent the accurate use of the solution for moderate or even reasonably small 
y (Y = 6 say). The large y results give two real roots, so ere obviously 
{ limited in usefulness. It is unlikely that 7 would be large enough to require 
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this solution. The large 7 solution is good down to y/16 = 3 or so, which 
is actually a very good range in the ] erturbation parameter; it just happens 
that for practical rotors Y falls far below this limit. Letting fi go to zero 
does not change the order of the aerodynamic coefficients (because of the 
constant terms) so the results for large and small Y should be good for 
all fi of order 1 or smaller. Letting fi go to infinity does change the order 
of terms in the analysis, so both the small Y and large Y results may be 


invalid for very large fi (above fi = 10 say). 


It would be very desirable to be able to use the small fi and small y 

results to construct composite root loci which are reasonably accurate for 

all values of fi and Y likely to be encountered in helicopter rotors. The 

small fi results would be used up to about fi = 0.5; then the small Y results 

would be used up to )* = 5 or so. The major obstacle to this is the lack 
o 

of an 0(7 ) analysis for the critical regions; the small Y analysis presented 
here, which was carried only to 0(7) in the critical regions, is not adequate 
for the accurate construction of loci for any except very small 7. 


Possible extensions of the solutions described here include the following: 

2 3 

a) Extend the small fi results to Ofyt ) for ImX =£, and to 0 (fi ) or 


Ofyi*) to handle the Im = 3/2 critical region. 


b) Extend the small 7 results to 0(7 ) in the critical regions. 


c) Extend the large y results to 0(y -1 ). 

d) Extend the large pi results to 0(pi~*). 

2 -1 

The most useful would be the 0(y ) and 0(pi ) results. A small y solution 
reliable to y/16 =^or so (which an 0(7 2 solution should accomplish) could 
be combined with the small jx solution to construct accurate composite 
loci. Together these results would then cover most of the range of pi and y 
of conventional rotors. The large pi case extended to 0 (pi, should be able 

to predict accurately the first instability boundary of the pi loci, which 
occurs at pi = 2 to 2. 5. These two cases are however also the ones involv ing 
the most work. 

APPLICATION OF PERTURBATION TECHNIQUES TO HELICOPTER DYNAMICS 

This section returns to the question of whether perturbation techniques 
might be profitably applied to more complicated or more realistic dynamic 
systems than the one considered here. As part of the answer, consider 
what these techniques will not do: obviously they can not give results for 
cases where there is no parameter that is either small or large, for example 
when y - 16 and pi = 1. However, the four cases considered together cover 
a good deal of the ranges of pi and y , and with primarily analytical results. 

For many helicopters the small p i case will be quite satisfactory alone. 

What the techniques can do also includes: 
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a) Since they give analytic solutions they provide more insight into 
the problem, as well as specific design criteria for the system; this feature 
is particularly important for nonlinear or time-varying systems, which 
have properties much different from those of constant coefficient, linear 
systems. 

b) Perturbation methods can find, and handle, cases that are very 
sensitive to the parameters, or that are difficult to solve accurately by 
numerical methods. 

c) The methods provide more insight into the rather unusual be- 
havior of the solution of periodic systems, by showing explicitly how the 
periodic coefficients modify the transient solutions and why they give tue 
root loci their characteristic behavior in the critical regions. 

d) Finally, even if the techniques are not used to find the complete 
solution, it only takes a little work to find out where the problems are 

(e. g. , critical regions and transition regions) and what the order of things 
is, which information would be of invaluable help in the numerical analysis 
of a system. 

The extension to more degrees of freedom or more realistic aero- 
dynamic coefficients would certainly make the analysis more complicated. 

In general however any study — analytic, computational, or experimental — 


T-V 
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of a system becomes more complicated as the accuracy of the modelling 
of the true system increases, and perturbation techniques are not expected 
to be an exception to this rule. Regardless of the system being studied, 
the position perturbation techniques occupy between simple linear analyses 
and complex nonlinear numerical calculations makes them a very powerful 
tool for providing both exact solutions and increased understanding of 
problems in rotor dynamics. 

The problems in rotor dynamics to which perturbation techniques 
might profitably be applied amount to all those involving nonlinear or 
periodic coefficients, and there are many of those. There is some additional 
work that might be done with the flapping dynamics problem (one degree of 
freedom), including for example 

a) teetering rotor; 

b) cantilever blade, with correct aerodynamic coefficients; 

c) Inclusion of stall and compressibility in the aerodynamics. 

2 

The solutions found in this paper might be extended to 0(7 ) for the small Y 
case and to 0(^1 for the large fi case. These extensions might prove 
very useful, or only slightly more so than the solutions to the order pre- 
sented here; but they should be examined for the simple single degree of 
freedom problem before being considered for more complicated systems. 




T 




Even the small /i solution could be taken a little farther — for example, to 
0( M 2 ) in the ImA critical region. Starting with two degrees of freedom, 
possible problems in coupled dynamics include 

a) flap dynamics of a gimballed rotor; 

b) pitch/flap dynamics; 

c) flap/lag dynamics. 

The pitch/ flap system is the problem of rotor flutter. The flap/lag system 
is particularly rich in possible variations; the problem has periodic co- 
efficients if n > 0 of course, but it f s also nonlinear (even in hover) due 
to the inertial coupling of the degrees of freedom. It is moreover very 
sensitive to blade root geometry, so that an articulated and a cantilever 
blade have quite different dynamic characteristics. Problems in coupled 
dynamics with more degrees of freedom include 

a) pitch/flap/lag dynamics (three degrees of freedom); 

b) Hap dynamics of an N-bladed rotor (N£3) with flapping feedback 
control in the fixed system (at least four degrees of freedom). 

While for these problems the advance ratio n would probably be of most 
value as a perturbation parameter, there will likely arise problems where 
other parameters are also useful. As long as a reasonable model is 
chosen for the system, and as much effort is given to the interpretation 
of the solution as to its derivation, perturbation techniques should prove 


quite useful in providing information about these problems, and many 
others in rotor dynamics and aerodynamics. 

This paper has demonstrated the methods of perturbation theory and 
has provided examples of the information about dynamic systems which 
may be obtained using them. The techniques have proved very useful 
for the problem studied. It should not be concluded however that the 
techniques presented are all there is to perturbation theory; there are many 
more methods that have not been touched on here. Perturbation theory is 
a powerful, and yet not very sophisticated, mathematical technique which 
should prove very useful in analyzing some of the problems of helicopter 


dynamics. 
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Appendix I. Eigenvalues of a Periodic System 

Consider a general system of differential equations with periodic 
coefficients; this may be reduced to a system of first order equations, and 
may therefore be written (in matrix notation) as 


-A A 
x = Ax 


where A(t) is periodic: A(t+T) = A(t). It may be shown that the solution to 
this differential equation can be obtained in the form 

x(t) =2^(0) eVujft) 

A A 

The are the eigenvalues; the eigenvectors are periodic, Uj(t + T) = u.(t); 
and the numbers q^(0) are constants obtained from the initial conditions. 

The theory that shows this is called Floquet theory. The solution in this 
form is a direct extension of the normal solution for a constant coefficient 
differential equation, which is charac*erized.by constant eigenvectors. 


The eigenvalues Xj may be obtained by the following procedure. 
The equation 




where 0(t) is a matrix, is integrated over one period, from t * 0 to t ® T, 
with initial conditions 0(0) = 1 (the unit matrix). Then if X c j are the 
eigenvalues of the matrix 

C » 0(T) 




th3 roots Xj are given by 


I 



or 



While the roots \ c j (as eigenvalues of a re al matrix C) must appear as real 
numbers or complex conjugate pairs, the eigenvalues are under no such 
restriction. The root loci of periodic systems are thus characterized by 
the type of behavior sketched below. 



-104- 


The horizontal portions of the loci can only appear at ImA = n/rev or 
n + \/ rev . If the parameter being varied, for example the advance ratio p , 
is such that at p = 0 the system is not periodic, then the roots at point A 
are complex corrugates. As p increases the periodicity of the system 
increases, and the roots start moving toward n/rev (or n + £ / rev) lines 
remaining complex conjugate pairs though) . At some critical p (the point 
P* on the locus) the loci reach ImA = n/rev, and then for still larger M 
the frequency remains fixed while the real part of one root is decreased 
and that of the other is increased This behavior should be compared with 
that of two roots of a constant coefficient system which start out as complex 
conjugates, meet at the real axis, and then proceed in opposite directions 
along the real axis. The existence of periodic coefficients in the equations 
of motion generalizes this behavior so that it can occur at any ImA = n/rev 

or n + $ /rev, not just ImA = 0. The property of the solution that allows this 
behavior is the fact that the eigenvalues Tij(t) are themselves periodic. 

For a single degree of freedom, second order system, let x R be the 
solution obtained from integrating the equation with initial conditions 
x(0) - 1, x(0) ■ 0; and let x p be the solution with initial conditions x(0) = 0, 
x(0) = 0. Then the roots are given by the quandratlc equation 

X C “ + c + *pCT) - 0 
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Appendlx TI. Solution of the Secular Equation 
The method of multiple time scales often leads to an ordinary 
differential equation of the form 

*||+(a + id)0 + (b + ic)0 = o 


where 0 is a complex quantity, and the constants a, b, c, and d are real, 
letting 


D 2 - d 2 - <b 2 + c 2 ) - d 2 - |(b - lo, | 2 


it may be verified that the solution of the above differential equation is 


D 2 X): 


0 _ e [A(d - D *■ i(b + irt)e 


ID# 


+ A(d + D + i(h + ic))e 


-iD^l 


where A is a complex constant 


D 2 - 0: 


0 - e” a ^ [A((d + i(b + ic))# + 1) 


+ B(d + i(b + ic»] 


where A arid B are real constants 


D 2 <0: 0 * e" a ^ [A(d + ID + i(b + ic))e D ^ 

+ B ;d - ID + i(b + ic))e -D ^) 
where A and B are real constants 
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The limiting case b - c = 0 gives D = d so the solution is 

0 = Ae -' a + td >* 
where A is a complex constant 

The region of decreased stability, i. e. , the region where the real part 

of the eigenvalue becomes more positive (the critical region), is given by 

2 2 
D <0. The boundary of the critical region is D - 0. One root in the critical 

region becomes less stable, but the other becomes more stable. Furlher- 

o 

more, inside the critical region (D“<0) there is a change in the real nart 

of the root but no change in the imaginary part, i.e. , the frequency; while 

2 

outside the region (D >0) there is a change due to D in the freque icy, but 

no change (hie to D in the real part. This behavior follows that expected 

2 

of the eigenvalues of periodic systems (see Appendix I). Indeed, D is 
a measure of the relative effects of the 0 and $ terns in the differential 
equation; the former usually comes from the constant coefficients in the 
equation of motion and the latter from die periodic coefficients. 
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Gradations on the loci are in tenths: 
except where noted 


/i= 0 .1 .2 .3 .4 .5 


ImX 







2 real roots, ImA s O 




Lines of constant ImX and ReX , based on the small n results 

(to order K p »0.1; ImX, R< 

circled values of ImX indicate areas in which ImX is constant. 
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Fig. 5. Lie 38 of constant ImX and ReX , based on the small M results 

o 

(to order y ); v = 1, Kp = -0. 1; —ImX , 

ReX , circled values of ImX indicate areas in which ImX is 
constant. 
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Fig ' 7 ‘ T ^ e averages of the 

01 “*• Mro ^ynaiBic coefficients, which give 

th© roots for small y 



(to order 7 ); v - 1 and K D = 1; the effect of the critical 
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Fig. 10. Motion of the root on the real axis for varying p , y 

For large y there are two real roots; x^ is the finite real 

root; the other root is at X ■ «• for Y * •». 
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Fig. 12. Sketch of the characteristic behavior of the critical region 

boundaries and stability boundaries for large p ; boundary 

of region in which Im is fixed at a/rev or n + — — 

boundary of region in which the real part of one root la 
positive, i. e. , unstable. 
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